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SUMMARY 


OBJECTIVE 

Analyze  techniques  and  technologies  for  performing  frequency  decomposition 
(i.e.,  Fourier  transformation)  on  an  input  signal  in  a  small-size,  low-cost,  low-power  unit. 
(Typical  device  parameters  of  interest  are  an  input  signal  bandwidth  of  2  MHz  and  a 
frequency  resolution  of  4  kHz.) 

RESULTS 

1 .  No  technology  satisfies  the  stated  objective  with  its  present  day  state-of-the-art. 

2.  Digital  FFT  circuits  using  high-speed  A/D  converter  and  multiplier  chips  meet 
the  small-size  and  low-cost  requirements  and  come  within  a  factor  of  4  of  the  speed  speci¬ 
fication  but  are  power  consumptive. 

3.  Incoherent  optics  technology,  particularly  the  electro-optical  processor  under 
development  at  NOSC,  and  electronic  CCD  technology,  particularly  the  chirp-Z  transform 
devices  developed  at  Reticon  and  Texas  Instruments,  both  have  high  potential  for  meeting 
the  stated  requirements  with  some  further  development. 

4.  The  remaining  technologies  investigated  are  clearly  inadequate  in  that  a 
major  technology  breakthrough  would  be  required  to  cost-effectively  meet  the  desired 
specification. 

RECOMMENDATIONS 

1.  Continue  monitoring  the  use  of  digital  electronic  technology  for  spectral 
analysis.  It  is  anticipated  that  increased  throughput  and  decreased  power  consumption 
should  be  achieved  in  the  near  future  due  to  the  considerable  momentum  of  the  commercial 
digital  electronics  field  and  the  advances  predicted  by  DoD’s  VHSIC  program. 

2.  Support  the  development  of  a  high-speed  DFT  module  using  electro-optical 
processing  (FOP)  technology.  The  use  ot  a  custom  gated-output  CCD  now  under  develop¬ 
ment  at  RCA  would  allow  the  fabrication  ot  a  128-point  DFT  analyzer  with  a  throughput 
rate  of  1 M  real  samples  per  second  in  FY  80  with  high  potential  for  further  performance 
improvement. 

3.  Support  the  development  of  high-speed  analog  circuitry  (either  on-chip  or 
external)  for  performing  the  pre-  and  post-multiplications  associated  with  a  chirp-Z  trans¬ 
form.  Electronic  CCD  devices  can  meet  the  desired  performance  specifications  for  the 
crosscorrelation  portion  of  the  algorithm  with  their  present  state  of  development. 
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1.  INTRODUCTION 


This  report  summarizes  the  results  of  a  project  entitled  “Multiple  Channel  Applica¬ 
tion  Study”  performed  by  the  author  for  Rome  Air  Development  Center/IRAA  during  the 
period  1  March  1978  through  30  September  1978.  The  following  four  paragraphs  constitute 
the  statement  of  work  for  this  task: 

OBJECTIVE:  To  investigate  and  report  on  the  feasibility  of  determining  activity  and 
demodulating  of  multichannel  communication  signals  (up  to  600  channels)  in  real  time. 

BACKGROUND:  In  recent  years  considerable  progress  has  been  made  in  technology 
and  the  use  of  this  technology  in  the  development  of  complex  communication  systems.  For 
high-capacity  systems  (greater  than  64  channels)  it  becomes  extremely  difficult  for  an  inter¬ 
cept  system  to  monitor,  track,  and  subsequently  demodulate  the  active  channels.  Recently, 
it  has  been  realized  that  there  may  be  important  potential  uses  for  optical  processing  techniques 
in  the  field  of  COMINT.  The  ability  to  perform  functions  such  as  filtering,  demodulation  and 
detection  at  extremely  high  rates  by  means  of  optical  techniques  could  make  possible  dramatic 
performance  improvements  and  in  some  cases  provide  entirely  new  capabilities  for  real-time 
signal  processing.  In  addition,  the  use  of  high-speed  solid  state  technology  offers  a  potential 
to  provide  solutions  to  the  multi-channel  communication  problem. 

EFFORT:  A  study  shall  be  conducted  to  determine  the  feasibility  of  real-time  pro¬ 
cessing  of  multi-channel  communication  signals  (up  to  600  channels).  Various  types  of 
technology  shall  be  investigated  to  determine  which  techniques  will  be  most  effective  in 
solving  the  problem.  These  technologies  shall  include  conventional  digital  electronics,  analog 
electronic  devices,  incoherent  electro-optical  processors,  and  optical  processors  utilizing 
fiber  optics.  Demonstration  of  the  selected  technique  or  techniques  is  desirable. 

DATA  ITEM:  A  report  detailing  the  technology  comparisons  and  recommending 
specific  approaches  shall  be  a  deliverable  requirement  of  this  effort. 

The  principal  signal  processing  problem  inherent  in  performing  the  desired  functions 
of  filtering,  demodulation,  and  detection  on  multichannel  communications  signals  is  the 
operation  of  spectral  decomposition.  Once  this  is  done,  the  channels  of  interest  can  each  be 
filtered,  demodulated,  and/or  detected  through  conventional  means.  The  more  mathematical 
term  for  spectral  decomposition  is  Fourier  analysis.  Hence  the  task  becomes  one  of  ( 1 )  inves¬ 
tigating  various  techniques  and  technologies  capable  of  performing  Fourier  transformations 
on  an  input  signal  in  a  small-size,  low-cost,  low-power  unit,  and  (2)  making  recommendations 
on  the  direction  of  further  technology  development  to  make  such  a  unit  feasible. 

To  clearly  illustrate  the  problem  limits,  the  following  typical  Fourier  transformation 
parameters  were  selected : 

-  an  input  signal  bandwidth  of  2  MHz, 

-  a  frequency  resolution  of  4  kHz  (i.e.,  the  input  signal  should  be  decomposed  into 
approximately  500  channels), 

-  the  use  of  a  good  window  function  (such  as  a  Kaiser-Bessel  window),  and 

-  an  appropriate  percentage  of  window  overlap  to  achieve  the  necessary  sampling 
rate  for  demodulation  (typically  50%). 

These  are  intended  simply  as  a  common  set  of  parameters  for  use  in  comparing  various 
techniques  and  states-of-the-art  and  are  not  meant  to  be  restrictive  or  necessarily  descriptive 
of  a  specific  application. 

The  next  section  will  define  the  Fourier  transform  operation  and  provide  a  mathe¬ 
matical  basis  into  which  we  can  insert  these  desired  parameters.  While  the  main  thrust  of 
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this  report  is  a  technology  comparison,  this  digression  into  theory  and  equations  is  necessary 
if  ambiguities  in  terminology  are  to  be  avoided.  With  that  task  completed,  the  subsequent 
sections  will  utilize  this  framework  in  quantitatively  assessing  the  ability  of  the  various 
signal  processing  technologies  currently  under  development  to  satisfy  these  demands. 

In  the  comparison  of  various  technologies,  several  specific  commercially  available 
devices  will  be  discussed  as  illustrations  of  the  state-of-the-art,  and  advertising  literature 
and  specifications  sheets  are  included  as  an  aid  to  the  reader.  This  should  not  be  interpreted 
as  an  endorsement  by  DoD  of  these  products  or  their  manufacturers. 

2.  THE  DISCRETE  FOURIER  TRANSFORM 

This  brief  section  on  the  Fourier  transform  utilizes  some  familiar  results  from  linear 
systems  theory  to  specify  a  set  of  parameters  to  be  utilized  in  subsequent  discussions.  No 
attempt  is  made  here  to  derive  any  of  these  results.  The  unfamiliar  reader  can  consult  any 
of  several  excellent  texts  on  linear  systems  theory,  digital  signal  processing,  and  Fourier 
analysis  (Refs.  1 ,  2,  3). 

The  continuous  Fourier  transform  is  a  relationship  between  a  function  of  time  (or 
space,  in  optical  systems)  and  a  function  of  temporal  frequency  (or  spatial  frequency). 

That  is,  given  the  input  function  g(t),  its  Fourier  transform  is 

G(f)  =  /°°  g(t)e“27riftdt  .  (1) 


Conversely,  given  G(f),  the  inverse  Fourier  transform  of  G(f)  is 
g(t)  =  J  G(f)  e27r 'Jtdf 


(2) 


The  unit  of  t  is  seconds  and  the  unit  of  f  is  cycles  per  second  (or  hertz).  These  equations  de¬ 
fine  a  relationship  between  two  different  ways  of  visualizing  a  signal:  either  its  conventional 
temporal  direct  waveform,  g(t),  or  its  temporal-frequency  representation,  G(f),  which  denotes 
the  amplitudes  and  phases  of  every  sinusoidal  component  of  g(t). 

In  more  graphical  terms,  Eqs.  (1)  and  (2)  relate  a  function  in  time  space  to  its  equiva¬ 
lent  function  in  frequency  space,  i.e., 


g(t) 


■t 


G(f) 


where  the  symbol  denotes  the  Fourier  transform  relationship. 

A.  COMPLEX  ANALYSIS 

Let  us  suppose  that  the  input  signal  to  the  receiver,  as  received  at  the  antenna,  is 
given  by  z(t);  i.e.,  it  is  a  real  continuous  temporal  waveform. 


The  Fourier  transform,  Z(t),  of  such  a  waveform  is  a  complex  continuous  waveform  in 
frequency  space.  It  is  Hermitian  (i.e.,  real  part  is  even,  imaginary  part  is  odd)  due  to  z(t) 
being  real.  Let  us  stipulate  that  due  to  some  a  priori  knowledge  of  the  nature  of  z(t),  only 
the  region  of  total  bandwidth  F,  centered  on  the  frequency  f  ,  is  of  interest. 


Z(f) 


Since  many  of  the  technologies  to  be  investigated  in  the  next  chapter  are  sampled-data  systems, 
the  input  signal  must  be  sampled  before  being  processed.  If  z(t)  is  sampled  directly,  then  the 
Nyquist  sampling  theorem  states  that  z(t)  must  be  sampled  at  at  least  twice  the  highest  fre¬ 
quency;  hence  the  sampling  frequency  would  have  to  be  at  least  2(fQ  +  F/2).  But  since  the 
total  bandwidth  of  interest  is  only  of  width  F,  this  is  obviously  inefficient.  It  makes  much 
better  sense  to  first  baseband  the  frequency  region  of  interest.  This  is  accomplished  by 
multiplying  z(t)  by  exp  (27rifQt),  which  has  the  effect  of  shifting  Z(f)  to  the  left  by  the  amount 
fQ.  Note  that  z(t)  exp  (27rif0t)  is  not  a  real  function  and  hence  its  Fourier  transform, 

Z(f-  fQ),  is  not  Hermitian. 


The  next  step  is  to  pass  the  basebanded  signal  through  a  low-pass  filter  that  transmits  only  the 
frequency  components  |  f  |  <F/2.  This  yields  a  waveform  which,  in  frequency  space,  looks 
like 


Sampling  this  waveform  at  the  Nyquist  limit  yields  a  sampling  frequency  of  2(F/2). 

The  term  “sampling  the  waveform”  means  multiplying  that  waveform  by  a  sequence 
of  Dirac  delta  functions  (normalized  by  the  width  At): 

OO 

At  5  (t  -  n  A  t)  . 

n=-oo 

By  the  convolution  theorem,  multiplication  in  time  space  implies  convolution  in  frequency 
space.  Therefore  the  Fourier  transform  of  the  waveform  must  be  convolved  with  the  Fourier 
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transform  of  the  delta  function  sequence,  which  can  be  shown  to  be  also  a  sequence  of  Dirac 
delta  functions,  namely, 


2  *('-£)• 

n=-oo 

Thus  sampling  this  time  waveform  at  intervals  of  At,  where 


has  the  effect  in  frequency  space  of  replicating  the  frequency  waveform  about  multiples  of  F. 


Also,  since  many  of  the  devices  to  be  discussed  have  sampled-data  outputs,  the  wave¬ 
form  in  frequency  space  needs  to  be  sampled.  Sampling  the  spectrum  at  intervals  of  Af  Hz/ 
sample  has  the  effect  in  time  space  of  restricting  the  input  signal  to  a  time  window  of  width  T 
seconds  and  replicating  this  function  every  T  seconds,  where 

Af  =  \  ■  (4) 

That  is 


Since  the  time-space  and  frequency-space  representations  both  completely  characterize  the 
same  physical  signal,  there  must  be  an  equal  number  of  “variables”  or  samples  in  each  repre¬ 
sentation.  One  can  think  ot  this  as  a  conservation  of  the  number  of  degrees  of  freedom.  Let 
this  number  be  N.  Thus 


_F 

Af 


(5) 


Since  by  sampling  in  time  and  frequency  the  processing  task  has  been  changed  from  a  con¬ 
tinuous  one  to  a  discrete  one,  the  integral  formulation  of  Eqs.  (1)  and  (2)  becomes  the  dis¬ 
crete  Fourier  transform  (DFT)  operation  of 
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and 


N 


where  n  and  m  denote  sample  numbers  in  time  space  and  frequency  space,  respectively. 

All  that  remains  of  this  analysis  now  is  to  insert  the  desired  typical  parameters  from 
Section  1  into  Eqs.  (3),  (4),  and  (5)  to  characterize  the  form  of  the  required  Fourier 
transformer. 

If  the  input  signal  bandwidth  of  interest  is 

F  =  2.048  MHz  (8) 

* 

then  by  Eq.  (3)  the  required  temporal  sampling  period  is 

At  =  --Sr  =  0.48828 1250 /as  •  (9) 

F 

If  the  frequency  resolution  of  interest  is 

Af  =  4  kHz  (10) 


then  by  Eq.  (5)  the  number  of  complex  samples  to  be  utilized  is 

_F  =  2.048  MHz 
N  Af  ,004  MHz 


(ID 


Inserting  these  two  results  into  Eq.  (4)  yields  the  total  integration  time  (or  window  length) 
required  as 

T  =  N  At  -  250 /as  .  (12) 


As  discussed  extensively  in  Ref.  4,  it  is  beneficial  to  apply  a  window  function,  w(t), 
to  the  input  signal  either  before  or  during  the  DFT  operation  in  order  to  reduce  the  frequency 
sidelobes  below  the  -13-dB  level  for  the  rectangular  window  (i.e.,  no  window  function  at  all) 
assumed  so  far  in  these  calculations.  Reference  4  lists  several  choices  and  quantitatively 
assesses  the  trade-offs  involved.  It  shows  many  plots  of  desirable  window  functions  w(t)  and 
their  Fourier  transforms  W(f)-  The  multiplication  of  the  input  data  by  one  of  these  window 
functions  effectively  gives  very  little  weighting  to  the  input  sample  points  near  the  ends  of 
the  window.  Hence,  it  is  necessary  to  use  some  overlap  factor  in  calculating  the  DFT’s  so  that 
all  input  samples  are  taken  into  account.  The  most  commonly  used  overlap  factor  is  50 7c. 

Using  the  parameters  just  calculated,  the  desired  processor  must  transform  a  250-/as  window  of 
the  input  signal,  gn,  every  1 25/us.  That  is,  every  125  /as  it  must  read  in  512  complex  input  samples 
and  generate  512  complex  output  spectral  samples.  The  effective  throughput  rate  is  therefore 

5 1 2  complex  samples  _  4  0%  ^  complex  samples  _  (l^ 

125  /as  second 

In  summary,  the  required  Fourier  transformer  has  as  its  input 

OO  OO 

gn  =  g(n  At)  =  z(t)  •  exp  (2n  ifQt)  •  w(t)  •  At  ^  5(t  -  n  At)  ©  5(t  -  mT)  (14) 

n=-oo  m=-°° 


7 


and  produces  as  its  output 


oo  oo 

Gm  =  G(m  Af)  =  Z  (f-fQ)  (g)  W(f)<g)  2  6(f-nF),A1'S  6(f"mAf>  <15> 

n=-oo  m=-°° 

where  the  symbol  (X)  denotes  convolution.  It  must  perform  this  operation  with  input  rates 
for  gn  and  output  rates  for  Gm  of  approximately  4  million  complex  samples  per  second. 

B.  REAL  ANALYSIS 


While  some  of  the  technologies  to  be  discussed  (e.g.,  coherent  optics  and  digital  elec¬ 
tronics)  are  capable  of  processing  the  complex  input  data  stream  described  above,  several 
others  (e.g.,  incoherent  optics  and  CCD’s)  are  restricted  to  real  input  signals  and  real  opera¬ 
tions.  Fortunately,  the  DFT  operation  can  still  be  performed  by  such  devices  by  processing 
the  two  real  components  comprising  each  complex  quantity  separately.  The  equations  can 
be  written  by  simply  inserting  the  identities 

exp  (-  ift)  =  cos  (2n  ft)  -  i  sin  (27t  ft)  (16) 

and 

exp(27rift)  =  cos  (2ff  ft)  +  i  sin  (2?r  ft)  (17) 


into  Eqs.  (1 )  and  (2)  and  requiring  that  the  input  signal,  g(t),  be  real.  The  latter  requirement 
can  be  insured  through  reconfiguring  the  receiver  processing  as  follows:  Instead  of  base¬ 
banding  the  received  signal  z(t)  as  before,  perform  the  following  three  operations: 

1.  Pass  z(t)  through  a  bandpass  filter  to  retain  only  those  frequency  components  of 
interest,  namely 


f0-3  «= 


f|<fo+i 


Multiply  this  waveform  by 


cos  2tt (f0 -f) 


3.  Pass  the  resulting  real  waveform  through  a  low-pass  filter  to  retain  only  those 
components  within 

-  F  <  |  f  |  <  +  F  . 


This  results  in  a  spectrum  of  bandwidth  2F  (twice  that  obtained  previously)  but  which  is 
Hermitian.  The  temporal  waveform  therefore  is  real. 


If  this  real  temporal  waveform  is  sampled  at  twice  the  highest  frequency,  that  is,  at 
intervals  of  At.  where 


_1_ 

At 


=  2F 


(18) 
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then  the  resulting  spectral  waveform  is  replicated  about  multiples  of  2F. 


Continuing  as  in  the  complex  analysis,  the  spectral  waveform  must  now  be  sampled  at  inter¬ 
vals  of  Af  Hz/sample.  As  before,  this  restricts  the  input  waveform  to  an  interval  of  T  seconds, 
where 


and  replicates  this  waveform  in  time  space  about  multiples  of  T  seconds.  That  is, 


COMPLEX 


The  conservation  of  degrees  of  freedom  argument  is  slightly  more  complicated  in  this 
analysis  than  before  because  of  the  need  to  consider  complex  variables  as  containing  two  real 
degrees  of  freedom.  Clearly  the  number  of  real  samples  in  the  time-space  window  is 

N'  (real)  =  -i-  (20) 

The  number  of  real  degrees  of  freedom  in  frequency  space  is  twice  the  number  of  complex 
samples;  thus 

N  (complex)  =  ^  (real)  (21) 


However,  the  Hermitian  constraint  implies  that  only  one-half  of  the  2F/Af  complex  samples 
are  independent;  therefore 


N  (complex) 


N'  (real)  =  J_  /2F\ 

2  2  y  Af  / 


Combining  Eqs.  (21 )  and  (22)  yields 


(22) 


N'  (real) 


J_  =  2F 
At  Af 


(23) 


This  differs  from  Fq.  (5)  by  the  factor  of  2  anticipated  by  the  conversion  to  real  data  from 
complex. 


MM 
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Proceeding  as  before,  if  the  input  signal  bandwidth  of  interest  is 

» 

F  =  2.048  MHz 


(8) 


then  by  Eq.  (18)  the  required  temporal  sampling  period  is 
At  =  =  0.244140625  ms  • 


If  the  frequency  resolution  of  interest  is 
Af  =  4  kHz 

then  by  Eq.  (23)  the  number  of  real  samples  to  be  utilized  in  the  DFT  operation  is 


N'  (real) 


2F  =  4.096  MHz 
Af  0.004  MHz 


1024 


(24) 


(10) 

(25) 


Inserting  these  two  results  into  Eq.  (20)  yields  the  total  integration  time  (or  window  length) 
required  as 

T  =  N' At  =  250  ms  •  (26) 

Again,  the  use  of  a  window  function  is  highly  desirable,  and  this  in  turn  necessitates 
the  use  of  some  overlap  factor  in  performing  these  DFT  operations.  A  50%  overlap  factor 
would  imply  that  the  desired  processor  must  transform  a  250-ms  window  of  the  input  signal, 
gn,  every  1 25  ms.  That  is,  every  1 25  ms,  it  must  read  in  1 024  real  input  samples  and  generate 
1 024  real  output  spectral  samples  (i.e.,  5 1 2  real  values  and  5 1 2  imaginary  values).  The 
effective  throughput  rate  is  therefore 

1024  real  samples  _  g  ^  real  samples  ^7) 

1 25  ms  second 

As  expected,  the  factor  of  2  difference  between  Eqs.  (13)  and  (27)  is  due  to  the  use  of  real 
instead  of  complex  input  samples.  The  processing  rate  (in  terms  of  real  operations  per 
second)  is  the  same  regardless  of  the  mathematical  notation. 

With  this  analysis  completed,  a  common  notation  has  been  developed  to  use  in  the 
following  technology  investigation. 


3.  FILTER-ARRAY  SPECTRUM  ANALYZERS 

In  seeking  ways  to  implement  the  desired  spectral  decomposition,  one  may  see  adver¬ 
tised  a  large  variety  of  “signal  analyzers,”  such  as  the  Hewlett  Packard  model  8565A,  described 
in  the  next  two  pages.  Such  units  perform  a  power  spectrum  analysis  on  incoming  signals 
with  frequency  content  of  several  MHz  to  several  GHz.  They  are  capable  of  achieving 
resolutions  as  Fine  as 

Af  =  1 00  Hz. 

They  are  readily  available,  utilize  existing  low-risk  technology,  are  compactly  packaged,  and 
arc  moderately  priced.  For  example,  the  model  advertised  in  the  next  two  pages  costs  about 
$1 8,000.00,  consumes  220  W,  and  occupies  IVi  in.  of  rack  space. 
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SIGNAL  ANALYZERS 

Spectrum  Analyzer,  10  MHz  to  40  GHz 

Modal  856SA 

•  0.01  to  22  GHz  coverage  with  internal  mixer 

•  Internal  preselection  1.7  to  22  GHz 

•  Wide  choice  of  resolution  handwidths 


•  Simple  three  knob  operation 

•  Absolute  amplitude  calibration 

•  CRT  bezel  readout  displays  control  settings 


8565A  Spectrum  Analyzer 

Covering  from  0.0 1  to  22  GHz  with  its  internal  mixer,  the  8565A 
has  built-in  preselection  and  brings  accuracy  plus  convenience  to  mi¬ 
crowave  spectrum  analysis.  The  wide  range,  spurious-free  display, 
compact  design  and  ease  of  use  make  it  well  suited  for  lab,  produc¬ 
tion,  or  field  applications  requiring  accurate  measurement  from  IF 
thru  microwave  frequencies.  The  8565A  can  cover  0.01  to  22  GHz  in 
just  two  spans  for  rapid  location  of  signals  prior  to  close-in  analysis  in 
one  of  six  bands.  Coverage  is  easily  extended  up  to  40  GHz  with  the 
HP  1 1 5 1 7A  External  Mixer. 

High  Resolution 

Fully  automatic  stabilization  in  narrow  spans  reduces  residual  FM 
and  drift.  Standard  resolution  bandwidths  range  from  I  kHz  to  3 
MHz.  The  1  and  3  MHz  bandwidths  allow  fast  sweeps  in  wide  spans 
and  increased  dynamic  range  for  pulsed  RF;  narrow  bandwidths  for 
measurement  of  closely  spaced  signals.  Option  100  provides  addition¬ 
al  100  Hz  and  300  Hz  IF  bandwidth  filters,  and  residual  FM  is  <100 
Hz  when  stabilized.  This  100  Hz  resolution  is  useable  up  to  8.5  GHz 
and  the  300  Hz  resolution  bandwidth  to  22  GHz.  All  resolution  filters 
are  gaussian-shaped  for  repeatable  measurements,  faster  non-distort¬ 
ing  sweeps  and  best  pulse  response. 

Absolute  Amplitude  Calibration 

Absolute  signal  levels  from  - 1 1 2  dBm  to  +30  dBm  are  easily  mea¬ 
sured  because  the  HP  8565A  always  displays  the  value  of  the  refer¬ 
ence  line  with  LED's  in  the  CRT  bezel  and  at  the  reference  level 
control.  Changes  in  RF.  IF  gain,  and  preselector  loss  are  automatical¬ 
ly  included.  In  addition,  flat  frequency  response  insjres  accuracy  for 
relative  as  well  as  absolute  power  measurements. 

Wide  Dynamic  Pange 

Internal  preselection  (1.7  to  22  GHz)  enables  you  to  measure  dis¬ 
tortion  products  as  much  as  100  dB  down.  Even  for  closely  spaced 
signals  or  measurements  below  1.7  GHz,  all  distortion  products  are 
greater  than  70  dB  down.  In  either  case,  maximum  dynamic  range  is 
assured  even  for  I  watt  signals  with  the  70  dB  input  attenuator.  An 
input  limiter  (0.01  to  1 .8  GHz)  and  the  internal  preselector  ( 1 .7  to  22 
GHz)  enable  the  8565A  to  withstand  RF  signals  up  to  +30  dBm  for 
all  input  attenuator  settings. 

Designed  For  Convenience 

Coupled  controls  allow  you  to  make  most  measurements  in  3  simple 
steps.  Green  color  coded  keys  preset  the  8565A  for  normal  operation 
so  a  measurement  only  requires  that  you  tunc  to  a  signal,  select  a 
desired  span,  and  raise  it  to  the  reference  level.  Automatically  select¬ 


ed  sweep  times  insure  a  calibrated  display  for  all  combinations  of  fre¬ 
quency  span,  resolution  bandwidth  and  video  filtering. 

The  CRT  bezel  LED’s  display  all  pertinent  control  settings  to  give 
you  all  the  information  needed  for  signal  evaluations  in  one  central 
location.  These  data  are  also  captured  in  CRT  photos. 

8444A  Option  H59  Tracking  Generator 

Make  swept  frequency  response  measurements  to  ±1.7  dB  from  1 0 
to  1 300  MHz  (  ±  2.5dB  up  to  1500  MHz)  with  greater  than  90  dB  of 
dynamic  range.  The  output  is  absolutely  calibrated  at  0  dBm  and  con¬ 
tinuously  variable  to  <-I0  dBm  .  The  frequency  of  unknown  signals 
as  well  as  the  frequency  of  any  point  on  the  frequency  response  curve 
can  be  measured  from  the  external  counter  output  using  the  low-c  »st 
HP  5300/5305B  Counter. 

8750A  Storage-Normalizer 

The  analyzer  is  made  even  easier  to  use  with  the  digital  storage  of 
the  8750A  because  there  is  no  need  to  re-adjust  intensity  or  persis¬ 
tence  as  the  sweep  time  changes.  With  the  push  of  a  button,  a  signal 
can  be  frozen  on  the  CRT  and  then  compared  directly  to  the  current 
input  signal.  Traces  can  also  be  compared  arithmetically  (i.e.,  nor¬ 
malized)  to  automatically  remove  frequency  response  variations.  This 
is  especially  useful  when  used  with  the  HP  8444A  Opt.  H59  Tracking 
Generator. 

8565A  Specifications 

Frequency  Specifications 

Frequency  range:  0.01  to  22  GHz  with  internal  mixer,  14.5  to  40 
GHz  with  HP  I I517A  External  Mixer. 

Tuning  accuracy  (digital  frequency  readout  in  any  span  mode) 
Internal  mixing:  0.01  to  2.5  GHz  <±5  MHz,  ±20%  of  Freq. 
Span/Div;  2.5  to  22  GHz.  ±0.2%,  ±20%  of  Freq  Span/Div. 
External  mixing:  14.5  to  40 GHz  <  ±0.7%,  ±  20%  of  Freq  Span/- 

Div. 

Frequency  spans 

1.7  to  22  GHi:  multiband  span  from  1 .7  to  22  GHz  in  one  sweep. 
Full  band:  displays  spectrum  of  entire  band  selected. 

Per  division:  eighteen  calibrated  spans  from  I  kHz  per  div.  to  500 
MHz  per  div.  in  a  I,  2,  5  sequence,  plus  a  full  band  span,  “F”. 

Span  width  accuracy:  the  frequency  error  for  any  two  points  on 
the  display  for  spans  from  500  MHz/div  to  20  kHz/div  (unstabi¬ 
lized)  is  less  than  ±  5%  of  the  indicated  separation;  for  stabilized 
spans  100  kHz/div  and  less,  the  error  is  less  than  ±  15%. 

Zero  span:  analyzer  becomes  a  manually  tuned  receiver. 
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Spectral  resolution  and  atabllity 

Raaolutton  bandwidth*:  resolution  (3  dB)  bandwidths  from  I 
kHz  to  3  MHz  in  1,3,  10  sequence.  Bandwidth  and  span  width  are 
independently  variable  or  may  be  coupled  for  optimum  display 
when  control  markers  are  aligned  (P4). 

Raaoiution  bandwidth  accuracy:  3  dB  points:  <  ±  15%. 
Salacttvlty  (60  dB/3  dB  bandwidth  ratio).  <15  1 
Stability:  Total  residual  FM  (fundamental  mixing  0.01  to  4.1 
GHz):  stabilized,  <200  Hz  p-p  in.  0.1  sec;  unstabilized  <10  kHz 
p-p  in  0.1  sec. 

Stabilization  rang#:  first  LO  automatically  stabilized  for  fre¬ 
quency  spans  100  kHz/div  or  less.  First  LO  residual  FM  typically 
30  Hz  p-p  when  stabilized. 

Noise  aldabanda:  >70  dB  down,  >30  kHz  from  center  of  CW 
signal  in  a  I  kHz  Res.  Bandwidth  and  a  10  Hz  (0.01 )  Video  Filter. 
Amplitude  Specifications 
Amplltuda  rang#  -  Internal  mixer 
Meaeurement  range 
Total  power:  +30  dBm  (1  watt). 

Damage  levels:  (501)  nominal  source  impedance.) 

dc:  0  V  with  0  dB  input  atten,  ±7  V  with  >  10  dB  input  atten. 

ac:  0  V  with  0  dB  input  atten.  10  V  peak  with  >10  dB  input 

atten. 

RF:  (signals  above  1 0  M Hz)  +  30  dBm  for  any  attenuator  setting. 
Gain  compreaalon:  <1  dB  for  0 dBm  input  level  with  OdB  attenu¬ 
ation. 

Average  noise  level:  max.  avg.  noise  level  with  I  kHz  Res.  Band¬ 
width  (0  dB  atten  and  3  Hz  video  filter)  is  in  the  table  below: 


Fraquancy 
Sand  (GHz) 

Flat  If 
in  MHz 

Harmonic 

Mods 

Horn  Laval 
(dBm) 

Fraquancy  Ratponta' 

( -  dB  MAX) 

18 

)050 

1- 

-11? 

12 

1741 

3)14 

1- 

-109 

17 

38-85 

3)14 

?- 

-103 

)5 

58-1)9 

3)14 

3- 

-94 

?5 

85-18 

3)14 

4+ 

-83 

30 

105-)? 

3)14 

5+ 

-75 

45 

‘Frequency  reeponee  include!  input  attenuetor.  preselector  and  mixer  frequency  response  plus  mix¬ 
ing  mode  gam  variation  (band  to  band) 


Amplitude  range  -  HP  11517A  External  Mixer 

Measurement  range  :  saturation  (gain  compression  <  I  dB).  - 1 5 
dBm.  Damage  level  >0  dBm  or  0.1  erg. 

Sensitivity  (Average  noise  level  in  a  10  kHz  IF  bandwidth): 
14.5-18  GHz  <-80  dBm,  18-26.5  GHz  <-70  dBm.  26.5-40 
GHz  <-60  dBm.  Typical  sensitivity  is  10  dB  better  for  each  band. 
Reference  Level 

Reference  level  range  +70  dBm  (  +  30  dBm  max.  input)  to  -102 
dBm  in  10  dB  steps  and  continuous  0  to -1 2  dB  calibrated  vernier. 
Reference  level  accuracy:  the  Auto  Sweep  setting  of  the  sweep 
time/div  control  insures  a  calibrated  display  within  these  limits: 
Calibrator  output  (100  MHz  ±  10  kHz):  -10  dBm  ±0.3  dB 
Reference  level  varl&llon  (Input  attenuator  at  O  dB).  10  dB 
steps  <  ±0.5  dB  (0  to  -70  dBm);  <  ±  1.0  dB  (0  to  -90  dBm). 

Vernier  (0  to  -12  dB)  continuous:  maximum  error  <0.5  dB. 
Input  attenuator:  (at  preselector  input,  0-70  dB  in  10  dB  steps) 
Step  size  variation:  <  ±  1.0  dB.  0.01  to  18  GHz;  <±1.5  dB, 
0.01  to  22  GHz. 

Maximum  cumulative  error  over  the  O  to  70  dB  range: 
<  ±  2.8  dB.  0.01  to  1 8  GHz;  <  ± 4.0  dB.  0.01  to  22  GHz. 

Frequency  response:  see  table  above. 

Switching  between  bandwidths:  3  MHz  to  I  kHz,  ±  1.0  dB 
Calibrated  display  range 

Log:  (  ’xpanded  from  reference  level  down):  70  dB  @  10  dB/div,  40 
dB  @  5  dB/div.  16  dB  @  2  dB/div  and  8  dB  @  1  dB/div. 

Linear:  full  scale  from  1 .8*i  V  (-102  dBm  in  501)  to  707  volts  (+70 
dBm)  in  10  dB  steps  and  continuous  0  to  -12  dB  vernier. 

Display  accuracy 

Log:  <±0.1  dB/dB,  but  <±1.5  dB  over  full  70  dB  display 
range. 

Linear:  <  ±01  division  over  full  8  division  deflection. 

Residual  responses  (no  signal  present  at  Input):  with  OdB  input 
atten.  fundamental  mixing  (0.01  to  4. 1  GHz)  <-90  dBm. 


Signal  identifier:  available  on  all  bands,  used  in  1  MHz/div  span  for 
signal  identification 

Signal  Input  Characteristics 
Input  501)  0.01  to  22  GHz 
Input  connector:  precision  Type  N  female. 

Input  Impedance 

Input  attenuator  at  0  dB:  50  ohms  nominal. 

SWR:  <1.5,0.01  to  1.8  GHz;  <2.0,  1.7  to  22  GHz  (at  analyzer 
tuned  frequency). 

Input  attenuator  at  10  dB  or  more:  50  ohms  nominal. 

SWR:  <1.3,  0.01  to  1.8  GHz;  <2.0,  1.7  to  22  GHz. 

LO  Emission  (2.00  to  4.46  GHz):  -50  dBm,  0.01  to  1 .8  GHz;  -85 
dBm.  1.7  to  22  GHz. 

Input  protection  (for  input  signals  from  0.01  to  22  GHz) 

0.01  to  1.8  GHz  frequency  band:  internal  diode  limiter. 

1.7  to  22  GHz  frequency  bands:  saturation  of  YIG  filter  (pre¬ 
selector)  occurs  at  total  input  signal  power  levels  below  input  mix¬ 
er  damage. 

External  mixer  input:  BNC  female  connector  is  a  port  for  LO  power 
transfer,  bias  current  and  IF  return. 

Sweep  Specifications 
Sweep  time 

Auto:  sweep  time  is  automatically  controlled  by  Frequency  Span/- 
Div,  Resolution  Bandwidth  and  Video  Filter  controls  to  maintain  an 
absolute  amplitude  calibrated  display. 

Calibrated  sweep  times:  21  internal  sweep  times  from  2  Msec/div 
to  10  sec/div  in  1,  2,  5,  10  sequence. 

Display  Characteristics 

Cathode  Ray  Tube  (aluminized  P31  phosphor,  8  x  10  dlv  inter¬ 
nal  graticule) 

Persistence 

Conventional:  natural  persistence  of  P3 1  phosphor. 

Write:  continuously  adjustable  from  0.2  sec  to  full  storage. 
Storage  time:  continuously  adjustable  from  I  minute  (full 
brightness)  to  >30  minutes  (minimum  brightness). 

CRT  Bezel  Readout:  bezel  LEDs  display  the  following  measurement 
data  (included  in  CRT  photographs  taken  with  the  HP  197A  Opt 
001, 006  Oscilloscope  Camera):  Ampl.  Scale  Factor,  Ref.  Level,  In¬ 
put  Atten.,  Res.  Bandwidth,  Sweeptime/Div.,  Freq.,  Freq.  Span/Div. 
General  Specifications 

Temperature  range:  operating  0°C  to  55°C,  storage  -40°  to 

+75°C. 

Humidity  range  (Operating):  <95%  R.H.  0°C  to  40°C 
EMI:  conducted  and  radiated  interference  is  within  the  requirements 
of  methods  CE03  and  RE02  of  MIL  STD  461A,  VDE  0871  and 
ClSPR  pub'n  1,  2  and  4. 

Power  requirements:  48-66  Hz,  100,  120,  200  or  240  volts  (-10% 
to  +5%)  220  V  A  max  (400  Hz  operation  available  as  Opt  400). 
Size:  188  H  x  426  W  x  552  mm  D  (7V**’  x  \6U”  x  21%"). 

Weight:  net  29.1  kg  (64  lbs).  Shipping  38.6  kg  (85  lbs). 

Standard  Options  Available 

Opt  100,  100  and  300  Hz  Resolution  Bandwidths:  adds  100  Hz 
and  300  Hz  resolution  bandwidths  with  11:1  shape  factor,  residual 
FM  <  100  Hz  when  stabilized  and  improves  sensitivity  by  10  dB. 

Opt  200— Calibration  in  dB^V 
Opt  400—400  Hz  Power  Supply 
Ordering  Information 
8565A  Spectrum  Analyzer 
Opt  100:  100  Hz  and  300  Hz  Resolution 
Bandwidths 

Opt  200:  Calibration  in  dB*iV 
Opt  400:  Internal  50  to  400  Hz  Power  Supply 
Opt  907:  Front  Handle  Kit 
Opt  908*  Rack  Flange  Kit 
Opt  909:  Rack  Flange  and  Front  Handle  Kit 
Opt  910:  Extra  Operating  and  Service  Manual 
11517A  External  Mixer  (taper  section  rcq’d) 

11518A  Taper  Section,  12.4  to  18  GHz 
115I9A  Taper  Section,  18  to  26.5  GHz 
11520A  Taper  Section,  26.5  to  40  GHz 
8444A  Opt  H59  Tracking  Generator,  10  to  1 500  MHz 
8750A  Storage-Normalizer 
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The  drawback  to  such  units,  at  least  for  the  application  of  interest,  is  that  they  com¬ 
pute  only  the  power  spectrum,  not  the  complex  amplitude  of  the  frequency  distribution. 

Hence  the  spectral  decomposition  cannot  be  used  for  complex  demodulation  as  required  herein. 


4.  DIGITAL  ELECTRONICS 
A.  GENERAL-PURPOSE  COMPUTERS 

In  a  typical  general-purpose  (GP)  computer  of  the  classic  Von  Neuman  architecture, 
input  data  is  brought  into  an  input  buffer  memory.  A  sequence  of  microprogram  instructions 
then  directs  the  arithmetic  logic  unit  (ALU)  to  perform  various  arithmetic  and/or  logic  opera¬ 
tions  on  this  data  and  store  it  back  in  memory.  When  the  desired  sequence  of  operations  is 
complete,  the  data  is  read  out.  A  software  program  to  directly  perform  the  multiplications 
and  summations  indicated  by  Eq.  (6)  can  readily  be  written;  thus  the  GP  computer  is  a  candidate 
to  be  considered. 

The  primary  differences  between  digital  computers  and  the  special-purpose  pipelined 
signal  processing  modules  to  be  discussed  shortly  (e.g.,  SAW  devices,  CCD’s,  optical  processors) 
are  in  versatility  and  speed.  Digital  electronic  components  can  perform  a  large  variety  of 
operations,  both  arithmetic  and  logic,  on  input  data.  (Hence,  the  primary  measure  of  perform¬ 
ance  is  “operations  per  second.”)  This  contrasts  sharply  with  the  analog  technologies  of 
interest,  which  can  only  perform  specific  operations,  such  as  multiplications.  As  we  shall  see, 
however,  the  increased  versatility  of  digital  computers  is  offset  by  much  lower  speed. 

The  DFT  algorithm  is  multiplication-intensive,  and  the  operation  of  “multiplication” 
is  one  that  requires  considerable  time  (compared,  say,  to  additions  or  shifts)  in  digital  imple¬ 
mentation.  Hence,  in  assessing  how  well  a  digital  electronic  processor  performs  the  DFT 
algorithm,  it  suffices  generally,  as  a  first  approximation,  to  consider  only  the  time  involved  in 
performing  the  required  multiplications  and  ignoring  the  processing  time  required  for  the  re¬ 
maining  operations. 

A  modern  state-of-the-art  minicomputer  system  can  perform  a  firmware  fixed-point 
multiply  instruction  in,  typically,  5  /us.  (For  example,  the  current  Hewlett-Packard  F-series 
CPU  performs  this  operation  in  6  /us.)  The  direct  computation  of  Eq.  (6)  requires  N^  complex 
multiplications  or  (N')2  real  multiplications.  The  discussion  of  Section  2B  shows  that  for  the 
parameters  of  interest,  the  window  length  is  1024  real  samples  and  the  required  throughput 
rate  is  one  1024-point  window  every  1  25  /us.  However,  a  simple  calculation  shows  that  a 
typical  general-purpose  computer  requires  about  5  s  to  perform  the  multiplications  for  each 
window,  i.e., 

d024)2  multiplications  t  5  /us  ^  5  s 

window  multiplication  window 

This  is  a  factor  of  40,000  too  slow. 

Since  the  limiting  factor  of  digital  electronic  implementations  for  signal  processing  tends 
to  be  the  restricted  multiplication  rate,  it  becomes  prudent  to  seek  alternative  algorithms  for 
calculating  the  Fourier  coefficients  [other  than  the  DFT  algorithm  of  Eq.  (6)]  that  reduce  the 
number  of  multiplications  at  the  expense  of  increasing  the  number  of  some  other  operation, 
which,  it  is  hoped,  the  digital  electronics  can  handle  quickly. 

Fortunately,  there  exists  such  an  algorithm  called  the  Fast  Fourier  Transform  (FFT), 
which  allows  the  .  mputation  of  the  spectral  values,  Gm,  in  only  N'  log->  N'  real  multiplications 
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instead  of  the  (N  r  real  multiplications  required  by  the  direct  implementation  of  Eq.  (6).  The 
price  paid  for  the  decreased  number  of  multiplications  is  a  complex  iterative  routing  of  data 
and  intermediate  results;  but  storing  and  retrieving  numbers  in  random-access  memory  is  a 
much  quicker  operation  than  multiplication,  so  the  use  of  this  algorithm  is  highly  beneficial. 
An  analysis  of  the  FFT  is  beyond  the  scope  of  this  report  and  is  well  described  in  Refs.  2,  3, 
and  5.  Using  the  FFT  algorithm,  a  GP  computer  with  a  5-#/s  multiplication  time  would  re¬ 
quire  50  ms  to  perform  the  multiplications  for  each  window,  i.e., 


1024  log2  1024 


multiplications 

window 


This  is  still  a  factor  of  400  too  slow. 


5  ms 

multiplication 


50  ms 
window 


Thus  software  or  firmware  implementation  of  the  spectral  analysis  operation  on  a  GP 
computer  is  far  too  slow  to  meet  the  desired  specifications. 


B.  ARRAY-TRANSFORM  PROCESSORS 


An  array-transform  processor  is  a  peripheral  unit  to  be  interfaced  with  a  GP  computer. 

It  is  specifically  designed  to  handle  those  operations  involving  large  arrays  of  numbers.  Most 
units  incorporate  a  hardware  multiplication  chip;  hence  they  are  considerably  faster  at  per¬ 
forming  the  FFT  operation  than  the  GP  computer  itself.  However,  they  still  retain  some 
degree  of  flexibility,  so  that  they  can  be  software-programmed  to  perform  other  array-inten¬ 
sive  and  multiplication-intensive  operations,  such  as  matrix  inversion  and  correlation  operations. 

As  typical  examples  of  the  state-of-the-art,  three  commercially  available  array  processors 
are  the  model  AP-1  20B  by  Floating  Point  Systems  Inc.,  the  model  100  by  Floating  Point 
Systems  Inc.,  and  the  model  MAP-300  by  CSP  Inc.  These  units  are  described  more  fully  in 
Refs.  6,  7,  and  8,  respectively.  The  only  specifications  of  interest  here  are  the  multiplication 
times,  which  are  167 , 250,  and  210  ns,  respectively.  Simple  calculation  shows  that  the 
fastest  unit  performs  the  multiplications  for  a  1024-point  window  of  real  input  data  in 
1 .7  ms,  e.g.. 


,  multiplications  167  ns  _  1.7  ms 

1 024  log-,  1  024  - r— -  •  - - - - - : ^  — — . 

z  window  multiplication  window 


(Actually,  the  AP-I  20B  specifications  sheet  lists  2.7  ms  as  the  time  required  for  a  1024-point 
real  FFT  when  the  time  required  for  all  the  other  operations  is  included.)  While  this  speed  is 
a  big  improvement  over  the  GP  computer,  it  is  still  a  factor  of  about  14  times  too  slow  for 
real  time  performance  of  the  task  stated  in  Section  2. 

The  accuracy  of  these  machines  is  very  good  due  to  the  use  of  floating  point  arithme¬ 
tic.  The  AP-1  20B  carries  38  bits  of  accuracy  (a  10-bit  binary  exponent  and  a  28-bit  two’s 
complement  mantissa),  while  the  MAP-300  uses  32-bit  numbers  (a  sign  bit,  a  7-bit  hexadecimal 
exponent,  and  a  24-bit  hexadecimal  mantissa). 


The  main  shortcoming  of  array  processors  for  the  intended  application,  however,  is 
not  their  speed  but  their  large  size,  large  power  consumption,  and  high  cost.  For  example,  the 
AP-1  20B  occupies  6.4  ft~\  consumes  1  200  W,  and  costs  approximately  S40K  (depending  on 
memory  options).  Additionally,  since  array  processors  cannot  operate  as  stand-alone  units, 
all  of  this  is  in  addition  to  the  size,  power  consumption,  and  cost  of  the  host  computer. 
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C.  SPECIAL-PURPOSE  DIGITAL  ELECTRONICS 

A  third  alternative  in  the  digital  electronics  arena  is  to  dispense  with  the  flexibility  and 
programmability  of  the  GP  computer  and  the  array  processor,  and  to  design  and  fabricate  a 
module  customized  to  the  FFT  algorithm  with  the  parameters  of  Section  2.  With  such  a 
special-purpose  system,  one  would  be  limited  only  by  the  speed  of  the  basic  multiplication 
operation.  (This  assumes  that  high-speed  bipolar  memory  is  used  so  that  the  memory-access 
time  is  not  the  critical  factor.)  Fortunately,  TRW  manufactures  a  line  of  very  high-speed 
multiplication  chips  well  suited  to  this  application.  Their  model  TDC  1 008  J  performs  the 
multiply-accumulate  operation  on  8-bit  inputs  in  70  ns,  and  their  model  TDC  1010  J  performs 
it  on  1 6-bit  inputs  in  1 1 5  ns.  Advertisements  for  these  two  products  are  shown  on  the  next 
two  pages  for  illustration. 

One  detail  of  the  FFT  algorithm  which  impacts  the  choice  between  these  two  chips 
is  that  round-off  errors  which  build  in  the  iterative  procedure  limit  the  output  accuracy.  A 
common  rule-of-thumb  is  that  one  loses  Vi  bit  of  accuracy  for  each  iteration.  Thus  a 
512-point  complex  FFT,  which  requires  log2  512  =  9  iterations  through  the  FFT’s  butterfly 
loop  loses  4 Vi  bits  of  accuracy.  If  8-bit  multiplications  are  performed,  the  final  answer  will 
have  only  8  -  4Vi  =  3Vi  bits  of  accuracy,  which  is  far  too  low  for  most  useful  applications.  Use 
of  the  16-bit  multiplication  chip,  on  the  other  hand  will  leave  16  -  4Vi  =  1 1  Vi  bits  of  accuracy. 
This  is  far  more  useful,  so  the  following  discussion  will  consider  this  chip. 

Assume  that  four  such  chips  are  utilized  and  are  configured  to  perform  a  complex 
butterfly  operation.  Then  the  time  required  to  perform  the  multiply-accumulate  operations 
for  performing  an  FFT  operation  on  a  512-point  window  of  complex  input  samples  is  0.53  ms, 
i.e., 


512  log2  512 


multiplications 

window 


1 15  ns  __  0.53  ms 
multiplication  ~  window 


This  is  only  a  factor  of  4  slower  than  the  1 25  ps  specified  in  Section  2A. 

An  estimate  of  the  size  and  power  requirements  of  such  a  system  can  be  obtained  from 
a  preliminary  circuit  design  for  a  stand-alone  FFT  module  performed  by  Dr.  Ed  Wrench  ot 
NOSC.  He  calculated  that  a  complete  system  would  require  about  75  chips  for  the  processing 
subsystem  and  about  25  chips  for  the  memory  subsystem.  These  100  chips  would  occupy 
two  5-  by  9-in.  circuit  cards.  The  4  multiplier  chips  would  require  3Vi  W  each;  the  14  high¬ 
speed  memory  chips  would  draw  ‘AW  each;  and  the  remaining  82  chips  would  average  about 
1/10  W  each.  Thus  the  total  power  consumption  would  be  about  30  W. 

The  speed,  small  size,  moderate  power  requirements,  and  relatively  good  accuracy  of 
this  approach  make  it  an  attractive  candidate  for  implementing  the  desired  operation. 
Furthermore  these  parameters  will  rapidly  get  better  due  to  (a)  the  considerable  momentum 
of  the  commercial  digital  electronics  field,  and  (b)  the  high  throughput  and  high  circuit  den¬ 
sity  promised  by  DOD’s  Very  High  Speed  Integrated  Circuits  (VHSIC)  program  in  the  signal 
processing  area. 


5.  INCOHERENT  OPTICAL  PROCESSING 

The  chief  advantages  of  optics  for  use  in  signal  processing  devices  are  (1 )  very  high 
multiplication  rates  and  (2)  inherent  parallelism  due  to  the  two-dimensional  nature  of  light 
wavefronts.  While  most  current  research  in  the  field  emphasizes  coherent  optics  (to  be  discussed 
in  the  next  section),  the  use  of  noncoherent  optics  has  much  to  offer.  The  price  paid  for 


y¥jwl  Now  you  can 
multiply.  •  •  or 
multiply  and  accumulate 

in70nsec 


TDC  1006J 

70  ns,  8  bits- $70  in  100’s 


■  Controllable  addition 
or  subtraction  in 
accumulator 

■  Round  control 

•  Bipolar  TTL  monolithic 
technology 

•  Power  dissipation  of 
1.2  watts 

•  Zero  hold  time 

•  Two's  complement  or 
unsigned  magnitude 


■  Accumulator  preloadable 
»  Cost  effective  as  a 

70  nsec  multiplier 

■  Multiply-accumulate 
in  70  nsec 

■  Ideal  for  complex 
multiplying  and  filters 
(including  FFTs) 

■  (Coming  soon:  16  bit 
multiplier/accumulator) 


Let  us  show  you  how  you  can  add  functions,  simplify  the 
design...  and  reduce  total  circuit  cost.  Available  from 
stock  from  Hamilton/Avnet  or  contact  your  local  TRW 
Electronic  Components  field  sales  office  or  call  us  at 
(213)  535-1831,  or  send  the  coupon. 
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TRW  takes  another  step 
in  digital  signal  processing 


Use  our  new  TDC  1010J  multiplier/ 
accumulator  (MAC)  to  build  a  high¬ 
speed  digital  signal  processor 
With  it  you  can  analyze  radar  signals 
or  X-ray  data;  communicate  with 
satellites  or  computers;  synthesize 
complex  waveforms — even  music. 

A  small  FFT  processor  based  on 
TRW's  new  MAC  operates  as  a 
spectrum  analyzer  too— add  one  to 
your  mini  or  micro  and 
you  don't^^^— 

sive  amounts  of  data  back  to  a 
number-crunching  mainframe  for 
reduction; you  can  reduce  it  right 
there  on  site  and  in  real  time ! 

It  can  analyze  voices,  earth¬ 
quakes,  geological  soundings  and 
submarine  signatures.  It  can  recog¬ 
nize  a  sticky  valve  in  an  automobile 
engine  or  in  a  human  heart 

There  was  a  time  when  the 
phrase  "FFT  processor"  conjured 
up  the  image  of  an  entire  bay  of 
sophisticated  electronic  hardware, 
but  that's  all  changed  now. 

Starting  with  just  a  single  TDC 


INTRODUCING 
A  16- BIT,  115-NSEC 
MULTIPLIER  /  ACCUMULATOR 


1010J,  you  can  design  your  own 
FFT  processor  on  a  small  pc  card 
It  will  operate  on  just  a  few  Watts  and 
the  CPU's  microcode  need  never 
even  touch  the  data. 

Simply  strobe  any  pair  of  16-bit 
numbers  into  the  MAC'S  on-chip 
input  registers  and  zip — the  chip 
delivers  the  correct  32-bit  product 
for  you  in  a  mere  115  nsec. 

An  on-chip,  35-bit  wide  accumu¬ 
lator  lets  you  choose  to  sum  a  series 
of  products  with  no  time  penalty  ; 
that's  both  a  double-precision  multi¬ 
ply  and  a  35-bit  add  in  the  same 
115  nsec  I 

Flexibility  is  a  key  feature  of  the 
TDC  1010J — it  works  on  numbers  as 
either  two's  complement  or  unsigned 
magnitude;  the  35-bit  accumulator 
can  be  directly  pre-loaded,  and  you 
can  round  off  the  accumulated 
products  to  single  precision. 

Our  new  MAC  is  fully  compatible 
with  industry  standards  TTL.flt  should 
be — after  all, TRW  invented  TTL 
and  patented  it  back  in  the  early  '60's, 
remember'’)  Of  course,  3-state 
output  buffers  are  provided 


TRW's  TDC  1010J  multiplier/ 
accumulator  is  packaged  in  a  64  pin 
DIR  It  consumes  just  3  Vi  Watts, 
uses  a  single  +  5V  supply  and  is 
radiation  hard.  It  is  priced  at  only 
$205  in  quantities  of  100 
Like  all  TRW  LSI  Products, 
the  TDC  1010J  multiplier/accumu¬ 
lator  is  available  now  from  stock 
through  Hamilton  Avnet  For  more 
information,  send  in  the  coupon  or 
talk  to  one  of  our  digital  signal 
processing  experts  at  213/535-1831 

TRW  LSI  Products  J 

An  Electronic  Components  Division  of  TRW  Inc 

P.0  Box  1125 

Redondo  Beach,  CA  90278 

Please  send  data  sheets  on  the  new  TDC  101 0J, 

16  x  16  bit  multiplier/accumulator 
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relinquishing  the  separate  control  of  amplitude  and  phase  information  obtainable  in  coherent 
systems  is  the  necessity  to  process  complex  data  as  two  real  quantities.  However,  as  discussed 
in  Section  2B,  this  does  not  present  any  difficulty  in  this  application.  What  one  gains  by 
utilizing  incoherent  light  is  the  freedom  from  vibration  problems,  thermally  produced  aberra¬ 
tions,  and  costly  high-quality  optics  characteristic  of  the  coherent  approach. 

The  two  papers  reprinted  in  Appendix  A  present  a  particular  approach  to  incoherent 
optical  processing,  named  the  electro-optical  processor  (EOP),  that  seems  well-suited  to  the 
problem  at  hand.  The  first  paper  (Ref.  9)  discusses  the  theory,  design,  and  application  of  such 
devices  and  ends  by  detailing  plans  for  fabricating  in  FY  80  a  CCD  that  would  give  an  EOP 
an  input  rate  capability  of  1M  samples/sec  while  performing  a  128-point  DFT  operation. 

The  second  paper  (Ref.  10)  gives  technical  details  of  presently  existing  EOP’s  capable 
of  performing  512-point  sliding-window  DFT  operations  with  input  rates  of  32K  samples 
per  second. 

The  strengths  of  this  approach  are: 

small  size  (2  in. 3  for  the  EOP  module  and  one  small  circuit  board  for  a  complete 
DFT  subsystem) 

-  low  power  (<5W) 

-  potentially  low  cost  (ultimately  <$1K  fora  512-point  1M  samples/sec  DFT  system) 

-  any  desired  window  function  can  be  incorporated  on  the  CCD  mask  and  processed 
with  no  extra  components  or  processing  time 

-  any  set  of  frequencies  can  be  calculated  (e.g.,  logarithmically  spaced  frequencies, 
or  only  over  a  limited  set  of  frequencies) 

-  useful  for  other  signal  processing  operations,  such  as  beamforming  and  correlation. 


6.  COHERENT  OPTICS  TECHNOLOGY 

Since  the  advent  of  the  laser  in  1960  the  field  of  coherent  optical  signal  processing 
has  attracted  considerable  attention.  Excellent  references  to  this  technology  are  the  text  by 
Goodman  ( Ref.  1  1 )  and  the  review  paper  by  Vander  Lugt  (Ref.  1 2).  As  mentioned  in  the  last 
section,  the  two  great  strengths  of  the  use  of  optics  (both  coherent  and  incoherent)  in  signal 
processing  are  (a)  very  high  multiplication  rates,  and  (b)  inherent  parallelism  due  to  the  two- 
dimensional  nature  of  light  wavefronts.  In  addition  to  these,  the  two  peculiar  strengths  of 
coherent  optics  are  (a)  the  complex  nature  of  the  light  field,  which  implies  that  the  amplitude 
and  phase  of  a  complex  number  can  be  propagated  through  the  system  as  the  amplitude  and 
phase  of  the  light  wavefront  at  a  given  point,  and  (b)  the  speed-of-light  calculation  of  the  two- 
dimensional  Fourier  transform.  In  a  coherent  optical  system  the  relationship  between  the 
complex  amplitude  light  distributions  in  the  front  and  back  focal  planes  of  a  lens  is  that  of  a 
Fourier  transform.  In  practice,  the  throughput  of  such  systems  is  limited  by  the  rate  at  which 
input  data  can  modulate  the  light  field  in  the  front  focal  plane  and,  to  a  lesser  degree,  by  the 
rate  at  which  the  output  values  can  be  detected  and  read  out  of  the  back  focal  plane.  Two 
types  of  real-time  input  transducers  of  interest  here  are  the  two-dimensional  spatial  light 
modulators  (SLM’s)  and  the  Bragg  cell.  These  shall  be  discussed  separately. 

A.  COHERENT  OPTICAL  PROCESSING  WITH  A  TWO-DIMENSIONAL  SLM 

Since  1960  considerable  effort  has  gone  into  the  development  of  real-time  two-dimen¬ 
sional  SLM’s  to  act  as  input  transducers  for  coherent  optical  systems.  The  desired  character¬ 
istics  of  such  devices  are: 
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-  read,  write,  erase  capability 

-  storage  capability 

-  high  frame  speed 

high  space-bandwidth  product 

-  grey  level  capability  (dynamic  range) 

-  high  uniformity 

-  long  recycle  lifetime 

-  environmental  stability 

-  economically  available. 

Several  candidate  materials  exist  and  are  currently  under  development.  These  include 

-  ferroelectric  ceramics 

-  liquid  crystals 

-  magnetic  bubble  domains 

-  Pockel’s  effect  crystals 

-  photodichroic  crystals 

-  deformation  devices 

-  membrane  light  modulators 

-  thermoplastics. 

Reference  1 3  by  Casasent  reviews  the  strengths  and  weaknesses  of  most  of  these. 

The  number  of  real  sample  points  per  DFT  window  in  the  specifications  of  Section  2B 
is  1024,  which  approximately  matches  the  number  of  resolution  cells  per  scan  line  of  these 
SLM’s  and  matches  the  one-dimensional  space-bandwidth  product  of  good  quality  optics.  So 
resolution  is  not  a  problem.  Also  the  speed  requirement  is  within  the  state-of-the-art.  Most 
of  these  SLM’s  can  operate  at  standard  television  rates,  namely,  63. S  ps  per  scan  line.  If  a  scan 
line  contains  the  1024  real  input  samples  then  this  technique  is  twice  as  fast  as  required  for  the 
desired  task. 

A  problem  with  coherent  optics  is  that  it  basically  is  not  a  good  match  for  the  appli¬ 
cation  under  discussion.  Coherent  optical  processing  systems  with  two-dimensional  SLM’s 
perform  two-dimensional  Fourier  transforms,  whereas  the  application  requires  only  one 
dimension.  True,  by  utilizing  cylindrical  optics  the  second  dimension  can  be  used  to  perform 
many  one-dimensional  transforms  in  parallel,  but  aspheric  optics  is  both  expensive  and  of  in¬ 
ferior  performance.  Also,  by  raster-scanning  the  input  and  output  planes,  the  second  dimen¬ 
sion  can  be  utilized  to  transform  a  very  long  input  signal,  but  the  1024  real  points  per  window 
of  interest  here  hardly  falls  into  that  category. 

A  more  serious  problem  with  this  approach  is  that  most  real-time  two-dimensional 
SLM  systems  are  bulky,  expensive,  and  not  ruggedized.  The  electron-beam  addressed  systems 
require  high  voltages,  and  the  laser-beam  addressed  units  require  registration  and  alignment. 
Major  breakthroughs  are  required  before  such  devices  become  practical  (i.e.,  rugged,  reliable, 
small  size,  low  cost,  low  power)  for  field  usage. 

B.  ACOUSTO-OPT1CAL  PROCESSING  USING  A  BRAGG  CELL 

As  depicted  in  Fig.  1 ,  a  Bragg  cell  is  a  device  for  converting  an  input  voltage  waveform 
into  a  bulk  acoustic  wave  propagating  through  a  crystal  material  that,  in  turn,  modulates  the 
phase  of  a  coherent  beam  of  light  passing  through  it.  It  is  a  very  useful  and  well-developed 
method  for  inputting  one-dimensional  real-valued  serial  data  into  a  coherent  optical  system. 
The  Applied  Technology  Divison  of  ITEK  Corporation  has  developed  several  such  systems  for 
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ACOUSTIC  ABSORBER 


Figure  1.  Acoustic  light  modulator. 

ECM  and  radar  signal  analysis  ( Ref.  1 4).  References  1 5  and  1 6  overview  some  recent  work  in 
the  field.  These  units  operate  with  center  frequencies  of  50  MHz  to  over  1  GHz  and  with  band- 
widths  of  from  a  few  MHz.  to  over  a  GHz.  Thus  the  total  bandwidth  is  more  than  required  for 
the  parameters  of  interest.  Also  the  1 000-spot  lime-bandwidth  product  matches  nicely  with 
the  1 024-point  windows  of  real  input  data  required.  The  problem  is  in  the  area  of  integration 
time.  Section  2B  requires  a  total  integration  time  of  250  /as  in  order  to  achieve  the  desired 
4  kHz  resolution.  Since  the  acoustic  velocity  in  lithium  niobate  is  6.5  mm//us,  a  Bragg  cell  of 
1.6  m  in  length  would  be  required.  Clearly  this  is  impractical.  Even  if  such  a  single  crystal 
could  be  grown,  the  acoustic  attenuation  from  one  end  to  the  other  would  be  a  problem,  and 
the  resultant  system  would  be  very  bulky. 

C.  COHERENT  DETECTION 

As  mentioned  earlier,  the  amplitude  and  phase  of  the  light  field  in  the  back  focal  plane 
of  the  lens  is  the  Fourier  transform  of  the  input;  however,  detecting  this  complex  light  field 
is  a  non-trivial  problem.  The  usual  system  simply  detects  the  intensity  (i.e.,  the  squared 
modulus  of  the  complex  light  values)  and  acts  therefore  as  a  power  spectrum  analyzer. 
Reference  17  discusses  techniques  in  which  a  reference  light  beam,  coherent  with  the  first,  is 
used  to  also  illuminate  the  output  plane.  The  intensity  of  the  resultant  interference  pattern 
then  yields  the  required  amplitude  and  phase  information.  Unfortunately,  optical  interfero¬ 
meters  are  notoriously  unstable,  hence,  extreme  precautions  must  be  taken  against  vibration, 
air  currents,  and  thermal  degradation. 


7.  FIBER  OPTIC  COMB  FILTER 

Another  signal  processing  technology  of  interest  here  is  the  use  of  fiber  optics  to  form 
a  transversal  filter.  The  basic  concepts,  as  reported  by  Taylor,  Dillard,  and  Hunt  (Refs.  18, 
19)  are  shown  in  Fig.  2. 
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FIBER 


Figure  2.  The  fiber  optic  transveral  filter. 

The  input  signal,  g(t),  temporally  modulates  the  optical  output  of  the  light  source, 
typically  an  injection  laser  or  a  light-emitting  diode  (LED).  (Since  the  intensity  out  of  an 
LED  is  linearly  related  to  the  current  in,  this  technique  is  capable  of  processing  analog  data.) 
This  light  is  then  coupled  into  a  bundle  of  N  optical  fibers  each  having  a  different  length,  Ln, 
where  n  =  1 ,  2,  3,  -  -  -,  N.  Thus,  each  fiber  introduces  a  different  delay  time  Tn  =  v  Ln,  where 
v  is  the  velocity  of  light  in  the  fiber.  As  shown  in  the  figure,  each  fiber  can  also  be  coupled 
through  an  optical  attenuator  in  order  to  modulate  the  intensity  transmittance  of  each  fiber, 
although  this  feature  will  not  be  utilized  here.  The  output  of  the  fiber  bundle  is  then  detected 
by  a  single  photodetector  where  a  summation  (i.e.,  energy  integration)  of  the  outputs  of  the 
N  fibers  takes  place.  It  can  be  shown  that  this  is  a  linear  system,  where  the  output  is  simply 
the  convolution  of  the  input  signal  with  an  impulse  response  defined  by  the  time  delays  and 
attenuation  factors  of  each  fiber.  In  other  words,  this  implementation  forms  a  finite  impulse 
response  (FIR)  filter  or  a  transversal  filter. 

If  each  fiber  delays  the  signal  by  some  amount  Tn  between  the  source  and  the  detec¬ 
tor  (i.e.,  no  attenuation),  and  if  the  length  of  each  fiber  is  an  integral  multiple  of  the  length 
of  the  shortest  one  (i.e.,  the  n^1  fiber  has  length  nL|),  then  the  impulse  response  of  the 
system  is  the  well-known  comb  function,  consisting  of  N  equally  spaced  impulses,  i.e., 

N 

h(t)  =  26(t~nT)  (28> 

n=l 


As  discussed  earlier,  convolution  in  time  space  is  equivalent  to  multiplication  in  frequency 
space  by  the  Fourier  transform  of  h(t).  The  Fourier  transform  of  a  sequence  of  N  equally 
spaced  impulses  is  a  series  of  sine  functions  of  mainlobe  width  47J-/NT  and  center-to-center 
spacing  lit II,  i.e., 
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OPERATE  IN  THIS  RANGE 


If  we  restrict  the  input  signal  to  contain  only  frequencies  between  7r/T  and  3n/T,  then  this 
“comb  filter”  acts  like  a  bandpass  filter  whose 

center  frequency  =  fc  =  ^  (29) 

and  whose 

bandwidth  =  BW  =  ^  .  (30) 

NT 


It  shall  be  mentioned  in  passing  without  further  elaboration  that  by  introducing  an  attenuation 
capability  into  each  fiber,  one  can  “shape”  the  impulse  response  to  obtain  filter  envelopes 
with  lower  sidelobe  levels  than  the  sine  function. 

The  velocity  of  light  in  fibers  is  nominally 


Thus  the  length  of  the  shortest  fiber,  the  one  producing  a  delay  of  T  seconds,  is  given  by 

Li  =  - - — —  meters,  (32) 

5  X  1 0“9 


and  the  length  of  the  longest  fiber,  the  one  producing  a  delay  of  NT  seconds,  is 


Liu  = - —  meters.  (33) 

5  X  1(T9 

Thus,  this  arrangement  of  N  fibers  constitutes  a  bandpass  filter.  One  can  easily  imagine  a 
second  group  of  N  fibers,  with  a  different  value  of  Lj,  leading  to  a  different  photodetector 
to  form  a  bandpass  filter  with  a  different  center  frequency.  It  is  a  straightforward  extension 
of  this  principle  to  construct  an  array  of  M  bandpass  filters  with  each  of  the  M  sets  of  N 
fibers  having  a  different  value  of  Lj.  If  the  center  frequencies  are  at  integer  multiples  of 
2tr/NT  across  the  range  of  7r/T  to  3n/T  Hz,  then  one  has  in  essence  a  spectrum  analyzer 
with  M  =  N  frequency  bins.  The  system  requires  one  LED,  N^  fibers,  and  N  photodetectors. 
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Let  us  now  consider  two  cases.  The  first  case  assumes  that  the  signal  z(t)  as  detected 
at  the  receiver  is  used  to  modulate  the  light  source  without  first  hederodyning  it  to  baseband. 
If  fc  is  1 09  Hz,  then 


T  =  71  =  seconds, 
c  109 

If  the  desired  bandwidth  is  BW  =  4  X  103  Hz  as  prescribed  in  Section  2,  then 


N  = 


27T 

(BW)T 


<4X,o3)(if) 


250,000. 


(34) 


(35) 


For  the  filtering  channel  centered  at  fc  -  2rr/T  Hz,  the  shortest  fiber  would  have  length 


h 


T  _  2n _ 

5  X  10~9  5  X  10'9  X  109 


1 .25  meters 


(36) 


and  the  250,000^  fiber  would  have  length 

l250,000  =  NL1  =  314.159  km.  (37) 

Clearly  such  a  large  number  of  fibers  (N2  =  6.25  X  1010)  of  such  long  lengths  (up  to  314  km) 
is  impractical. 

Case  number  two  assumes  that  the  received  signal  is  heterodyned  down  to  a  center 
frequency  of  2  MHz  before  being  inserted  into  this  transversal  filter.  If  fc  is  2  X  10^  Hz,  then 

T  =  -t—  =  - —  —  seconds.  (38) 

fc  2  X  106 


If  the  desired  bandwidth  is  still  4  kHz,  then 


N  = 


2tr 


2n 


(BW)T 


«4xio3»(2xToo) 


=  500. 


(39) 


(fc  was  chosen  in  order  to  make  N,  the  number  of  frequency  bins  observable,  come  out  to 
approximately  5 1 2  as  desired  for  the  parameters  of  Section  2.)  For  the  filtering  channel 
centered  at  f  =  27r/T  Hz,  the  shortest  fiber  would  have  length 

Li  =  - - —  =  - - tt  ~  628  meters  (40) 

5  X  10-9  (2  X  106)(5  X  10“9) 

and  the  500^  fiber  would  have  length 

l500  =  NL1  =  314  km-  (4,) 

Thus,  for  a  spectrum  analyzer  system  to  meet  the  specifications  desired,  one  would  require 
N2  =  250,000  fibers  ranging  in  length  up  to  314  km.  It  is  the  author’s  opinion  that  while 
such  a  unit  could  be  fabricated  for  laboratory  testing,  major  breakthroughs  in  the  art  of  fiber 
bundle  packaging  are  required  to  make  such  an  implementation  size  and  cost  effective. 
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While  this  sytem  has  been  discussed  in  terms  of  power  spectrum  analysis,  it  can  be 
extended  to  perform  a  complex  Fourier  transform  either  by  (a)  using  four  such  units  in 
parallel  (two  each  for  the  in-phase  and  quadrature  inputs)  with  the  impulse  response  of  each 
group  of  N  fibers  tailored  to  be  sine  or  cosine  waveforms  of  various  frequencies,  or  (b)  using 
phase  and  amplitude  modulation  of  a  carrier  as  the  input  and  complex  detection  to  remove 
the  effects  of  the  carrier  at  the  output.  Such  modifications  obviously  add  to  system  complexity. 


8.  IMPLEMENTATIONS  OF  THE  CHIRP-Z  TRANSFORM 


Equation  (6)  is  the  traditional  form  of  the  discrete  Fourier  transform  equation.  How¬ 
ever,  if  the  following  substitution  is  made: 


e~2mn  =  exp  [m^  -  (n  +  m)2  +  n^]  , 


(42) 


one  arrives  at  an  alternative  form  of  the  Fourier  transform  known  as  the  chirp-Z  transform, 
CZT,  (Ref.  3),  namely 

N 


Gm  "  cx'’ 


(»tt)  ±<M4) 

n=l 


exp 


■  (n  +  m)^ 
N 


(43) 


This  suggests  the  following  algorithm  for  performing  the  Fourier  transform  operation: 

1.  Multiply  ffi  by  the  chirp  waveform  exp 

2.  Cross-correlate  this  result  with  the  chirp  waveform  exp  for  a  number  of 

shifts  m,  and 


3.  Multiply  this  result  by  the  chirp  waveform  exp  ( i7r  -!2- 

\  N 

The  resulting  coefficients  Gm  obtained  from  the  CZT  algorithm  are  identical  to  those  obtained 
using  the  traditional  algorithm  of  Eq.  (6)  even  though  quite  a  different  mathematical  procedure 
is  used. 

The  use  of  this  algorithm  can  be  very  beneficial  if  used  in  conjunction  with  hardware 
that  performs  the  cross-correlation  operation  at  very  rapid  rates.  Two  such  hardware  imple¬ 
mentations  arc  (a)  surface  acoustic  wave  (SAW)  devices,  and  (b)  charge-coupled  devices 
(CCD’s). 

A.  SURFACE  ACOUSTIC  WAVE  DEVICES 


SAW  devices  represent  a  very  attractive  implementation  for  the  cross-correlation  por¬ 
tion  of  the  CZT  operation.  References  20  and  21  give  more  information  than  can  be  presented 
here.  These  devices  work  with  either  analog  or  sampled  analog  inputs  and  outputs.  A  typical 
device  consists  of  a  substrate  of  piezoelectric  material  on  one  face  of  which  is  deposited,  by 
photolithographic  techniques,  a  pattern  of  aluminum  or  other  conductor  electrodes.  In  its 
simplest  configuration,  these  sets  of  interdigitated  finger  electrodes  are  spaced  at  the  sampling 
rate  distance  through  the  use  of  the  relationship  ds  =  CRts,  where  ds  is  the  tap  spacing,  CR  is 
the  Rayleigh  wave  velocity,  and  ts  is  the  sampling  increment.  For  a  typical  substrate  of 
ST-cut  quartz,  CR  3  mm/gs.  Thus  for  an  integration  time  of  250  /as  (necessary  to  provide 
a  resolution  ot  4  kHz),  a  quartz  crystal  of  length  0.75  m  would  be  required.  This  is  far  be¬ 
yond  the  current  capability  ( i.e. ,  several  inches):  hence  SAW  technology  is  not  a  good 
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candidate  to  meet  the  specifications  of  Section  2.  For  applications  requiring  sampling  rates 
of  from  10  to  100  MHz,  with  carrier  frequencies  of  from  50  to  500  MHz,  and  with  typical 
fractional  bandwidths  of  about  10%  (i.e.,  5  to  50  MHz),  SAW  transversal  filters  for  implemen¬ 
tation  of  the  CZT  hold  great  promise  for  low-cost,  small-size,  low-power  systems. 

B.  CHARGE  COUPLED  DEVICES 

A  second  technology  worth  considering  for  implementation  of  the  cross-correlation 
portion  of  the  CZT  algorithm  is  charge  coupled  devices  (CCD’s).  In  fact,  the  advertisement 
copied  on  the  following  page  shows  a  commercially  available  CCD  chip  designed  for  just  that 
purpose.  References  22,  23,  and  24  give  more  information  than  can  be  presented  here. 

Basically,  CCD’s  are  sampled  data  analog  circuits  which  can  be  fabricated  by  metal 
oxide  semiconductor  (MOS)  technology  as  LSI  components.  As  such  they  are  directly  com¬ 
patible  with  other  MOS  circuits.  State-of-the-art  CCD  transversal  filters  have  operated  as 
video  devices  with  sample  rates  up  to  10  MHz.  CCD’s  operate  by  the  manipulation  of  injected 
minority  carriers  in  potential  wells  under  MOS  capacitors  and  thus  behave  as  capacitive 
reactances  with  low  power  dissipation.  However,  the  potential  wells  that  contain  the  minority 
carriers  also  attract  thermally  generated  minority  carriers;  hence  there  is  a  maximum  storage 
time  for  the  analog  signals  that  depends  on  the  dark  current  associated  with  the  temperature 
of  the  silicon.  At  room  temperature,  storage  times  of  hundreds  of  milliseconds  are  typical. 
Several  different  techniques  are  used  to  achieve  unidirectional  charge  transfer.  The  first  de¬ 
veloped  was  a  three-phase  clocking  structure  in  which  the  direction  of  charge  propagation  is 
determined  by  the  sequence  of  potentials  applied  to  the  three  electrodes. 

In  the  so-called  surface  channel  CCD’s,  if  the  minority  carriers  are  allowed  to  collect 
at  the  semiconductor-oxide  boundary,  poor  charge  transfer  efficiency  will  result  due  to 
minority  carriers  getting  caught  in  trapping  sites.  This  means  that  the  CCD  will  behave  non- 
li nearly  unless  there  is  sufficient  propagating  charge  present  to  fill  all  of  the  traps.  By  biasing 
the  operating  condition  of  the  CCD  so  that  about  10%  of  the  dynamic  range  is  used  for  the 
injection  of  a  “fat  zero”  bias  charge,  the  traps  are  kept  continuously  filled  and  the  device  has 
over  a  60  dB  dynamic  range.  In  practice,  the  input  signal  to  be  processed  is  added  to  a  fixed 
bias  somewhat  larger  than  one-half  of  the  peak-to-peak  value  of  the  signal  before  it  is  inserted 
into  the  CCD. 

In  general,  CCD’s  have  throughput  rates  of  from  10  kHz  to  10  MHz.  For  a  500-point 
line-array  CCD,  this  translates  to  an  integration  time  window  of  from  50  ms  to  50  ps,  which, 
for  a  Fourier  transform  operation,  translates  to  frequency  bin  resolutions  of  20  Hz  to  20  kHz. 
These  throughput  rate  and  resolution  figures  indicate  that  CCD  technology  is  a  good  candidate 
for  meeting  the  desired  specifications. 

So  far,  we  have  only  discussed  the  cross-correlation  portion  of  the  CZT  algorithm, 
with  the  result  that  CCD  technology  shows  considerable  promise.  We  have  yet  to  discuss  the 
pre-  and  post-multiplication  portions  of  the  CZT  operation  and  the  generation  of  the  required 
chirps  waveforms.  One  approach  to  these  tasks  is  to  use  off-the-shelf  digital  electronic  technol¬ 
ogy,  namely: 

1.  Perform  an  A/D  conversion  on  the  incoming  signal. 

2.  Generate  digital  versions  of  the  chirp  waveforms  by  reading  out  high-speed  ROM’s. 

3.  Use  the  aforementioned  TRW  high-speed  multiplication  chip. 

4.  Perform  a  D/A  conversion  to  obtain  an  analog  signal  to  insert  into  the  CCD 
cross-correlator. 


The  spectrum 
analyzer  on  a  chip. 

Where  else  but  from  Reticon. 

What  was  previously  thought  impossible  will  now 
be  an  everyday  occurence.  Moving  up  to  the  next 
stage  of  complexity  in  CCD  devices  has  resulted  in 
the  Reticon  R5601,  a  512  point  Discrete  Fourier 
Transformer.  This  technology  offers  a  spectrum 
analyzer  with  small  size,  light  weight,  low  power, 
high  reliability,  and  a  remarkable  low  cost.  Along 
with  its  associated  circuitry,  it  performs  the  Chirp 
Z  algorithm  to  give  a  256  spectral  line  display  in 
less  than  250  /-tsec.  It’s  small  enough  to  fit  into 
your  system,  yet  powerful  enough  to  have  a  signal- 
to-noise  ratio  in  excess  of  70db.  The  numerous 
applications  possible  include  speech  recognition, 
target  identification,  vibration  analysis,  bandwidth 
compression,  communications,  and  general  signal 
analysis. 

Currently  available  is  a  self-contained  evaluation 
module  on  two  printed  circuit  cards  just  80  square 
inches.  Just  hook  up  your  ±20  volts,  display  and 
you’re  on  the  air.  Use  the  on-board  oscillator  or 
externally  control  the  sampling  rate. 

The  R5601  is  the  latest  in  our  growing  family  of 
discrete  time  analog  signal  processing  devices.  All 
available  through  our  worldwide  network  of  over  20 
distributors  and  more  than  70  salesmen. 

Discover  the  1C  that  does  it  all. 

RETICON* 

AN  COMPANY 

910  Benicia  Ave  •  Sunnyvale,  California  94086 
(408)  738-4266 
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5.  Perform  an  A/D  conversion. 

6.  Repeat  steps  2  and  3  for  the  post-multiplication  operation. 

7.  Perform  a  D/A  conversion  to  provide  an  analog  output  for  display. 

This  approach  would  consume  too  much  power  to  really  be  cost-effective.  An  alternative 
approach,  which  eliminates  the  A/D  and  D/A  converters,  is  the  use  of  multiplying  DAC’s, 
wherein  a  digital  signal  (i.e.,  the  chirp  waveforms  from  ROM)  controls  the  gain  of  the  analog 
input  signal.  Unfortunately  state-of-the-art  multiplying  DAC’s  have  maximum  throughput 
rates  of  a  few  hundred  kHz,  which  is  slower  than  the  8M  samples/second  rate  specified  in 
Section  2. 

The  use  of  analog  technologies  to  perform  these  operations  in  a  cost-effective  manner 
appears  promising.  Texas  Instruments  has  demonstrated  a  32-point  C'ZT  device  with  the 
chirp-generators,  multipliers,  and  CCD  cross-correlator  all  on  the  same  silicon  chip.  An  alter¬ 
nate  approach  would  be  the  development  of  high-speed  switch-capacitor  multiplying  D/A 
converter  chips,  and  the  generation  of  analog  chirp  waveforms  by  pulsing  a  CCD  transversal 
filter. 

The  CCD  CZT  approach  has  promise  and  appears  a  fruitful  area  for  further  develop¬ 
ment.  Let  us  now  look  at  what  is  commercially  available  today. 

Reticon  Corp.  (Sunnyvale,  CA)  currently  sells  a  single  circuit  board  CZT  device 
employing  a  CCD  convolver  chip  for  about  $1600.  It  performs  only  a  256-point  power  spec¬ 
trum,  has  a  throughput  rate  of  only  100  kHz,  a  dynamic  range  of  70  dB,  an  accuracy  of  only 
8%,  and  computes  each  frequency  from  a  different  sliding  window  of  the  input  data.  Thus, 
their  present  device  does  not  meet  the  specifications  of  Section  2. 

This  circuit  board  could  be  straightforwardly  modified  to  perform  a  512-point 
Fourier  transform  (i.e.,  real  and  imaginary  inputs  and  outputs)  in  addition  to  perform¬ 
ing  the  power  spectrum,  have  a  throughput  rate  of  200  kHz,  and  an  accuracy  of  1%.  The 
speed  is  limited  by  the  analog  circuitry  for  the  pre-  and  post-multiplications.  The  throughput 
rate  could  be  increased  to  2  MHz  by  employing  digital  circuitry  at  the  expense  of  increased 
power  consumption  and  circuit  complexity,  as  discussed  earlier.  Also,  the  sliding-window 
feature  could  be  circumvented  by  employing  a  recirculating  delay  line  in  the  input  at  the 
price  of  halving  the  throughput  rate  or  halving  the  number  of  points  in  the  transform.  Alter¬ 
natively,  the  throughput  rate  could  be  increased  by  the  development  of  high-speed  analog 
multiplying  chips  so  that  a  throughput  rate  of  about  1  MHz  could  be  achieved  without  the 
power  and  expense  of  digital  external  circuitry. 


9.  CONCLUSIONS 

No  technology  meets  the  specifications  of  Section  2  with  its  present  day  state-of-the- 
art.  Both  the  electro-optical  processor  (EOF)  technology  (performing  a  DFT  or  FIR  filter  bank) 
and  the  electronic  CCD  technology  (performing  a  CZT)  have  high  potential  for  doing  so  with 
some  further  development.  The  EOP  requires  the  development  of  a  512  by  512  resolution-cell 
CCD  with  a  shift-register  rate  of  4  MHz.  Four  such  devices  would  then  be  required  to  meet  the 
specifications  of  Section  2.  Such  developments  are  currently  in  progress  at  the  Naval  Ocean 
Systems  Center  and  at  RCA  Laboratories  (David  Sarnoff  Research  Center)  funded  by  the 
Naval  Electronics  Systems  Command  (NAVELEX).  The  electronic  CCD  CZT  processor  requires 
the  development  of  high-speed  analog  multiplier  technology  (either  on  or  off  the  CCD  convol¬ 
ver  chip).  These  efforts  are  currently  being  pursued  at  both  Reticon  and  Texas  Instruments. 
Due  to  the  promising  outlook  offered  by  the  electro-optical  processor  and  electronic  CCD 
technologies,  more  information  on  them  is  provided  in  Appendices  A  and  B,  respectively. 
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Digital  FFT  circuits  employing  TRW’s  high-speed  A/D  converter  and  multiplier  chips 
meet  the  small-size  and  low-cost  requirements  and  come  within  a  factor  of  4  of  the  speed 
specification  but  are  power  consumptive  (say  30  W  compared  to  just  a  few  watts  for  the  EOP 
and  CCD  devices).  It  is  anticipated  that  increased  speed  and  decreased  power  consumption 
should  be  achieved  in  the  near  future  due  to  the  considerable  momentum  of  the  commercial 
digital  electronics  field  and  the  advances  predicted  by  DOD’s  VHS1C  program. 

The  remaining  technologies  are  clearly  inadequate,  i.e.,  a  major  technology  break¬ 
through  would  be  required  to  cost-effectively  meet  the  desired  specifications. 
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Introduction 

The  authors  have  developed  a  new  electro-optical  signal 
processing  module,  just  a  lew  cubic  inches  in  size,  capable  of 
performing  a  broad  variety  of  general  linear  filtering  opera¬ 
tions  at  very  rapid  rates  -  typically  greater  than  I0‘*  multipli¬ 
cations  per  second.  The  distinct  feature  of  this  processor  is 
its  use  of  incoherent  optics  to  achieve  speed  and  parallelism 
in  its  operation.  References  1-6  describe  previous  work  in 
this  technology.  The  present  paper  reviews  the  characteris¬ 
tics  of  this  processing  technique,  describes  the  mathematical 
operations  it  is  capable  of  performing,  gives  the  operating 
parameters  of  the  current  implementation,  and  presents  plans 
for  modifications  to  even  further  improve  its  performance. 

Processor  Description 

A  mathematically  rigorous  treatment  of  the  concept 
underlying  the  processor  under  discussion  is  contained  in 
references  4  and  5  :  however,  for  overview  purposes,  the  fol¬ 
lowing  qualitative  description  will  be  presented  here.  The 
processor,  illustrated  in  Figure  I.  consists  of  three  compo¬ 
nents:  a  light-emitting  diode  (LED)  illuminating  a  photo¬ 
graphic  transparency  (or  matrix  mask)  which  in  turn  is  in 
contact  with  an  area-array  charge-coupled  device  (CCD). 

Let  the  input  signal  be  an  analog  sampled-data  sequence, 
denoted  by  fn0+n  with  a  sample  rate  of  Vjn  samples/sec. 

Here.  nQ  is  used  to  denote  an  arbitrary  discrete  time  origin. 
This  signal  is  used  to  temporally  modulate  the  optical  output 


Figure  I.  Basic  components  of  the  electro-optical 
signal  processor. 


of  the  LED.  Thus  the  irradiance  illuminating  the  photo¬ 
graphic  transparency  is  directly  proportional  to  the  input 
signal  'n0+n  ant*  's  spatially  uniform  across  the  entire  surface 
of  the  mask.  The  transparency  contains  an  M  by  K  array  of 
elements  with  optical  transmittance  values  hm y  arranged  in 
a  geometrical  pattern  and  scale  identical  to  that  of  the  photo¬ 
site  array  in  the  CCD  *  Since  the  irradiance  transmitted  by  a 
transparency  is  the  product  of  the  illuminating  irradiance 
times  the  intensity  transmittance  of  the  mask,  the  light  dis¬ 
tribution  immediately  behind  the  m.k1*'  element  of  the  trans¬ 
parency  at  the  time  when  the  n0+n'*1  input  sample  is  modu¬ 
lating  the  LED  is  proportional  to  fn  +„  bn,  ^  Assuming 
that  the  mask  and  the  optically  sensitive  area  of  the  CCD  are 
in  physical  contact  and  registered  such  that  the  I’.i.k1 '  cell  of 
the  transparency  is  superimposed  on  the  nt.k1*'  element  of 
the  CCD.  then  this  irradiance  pattern  is  detected  by  the  CCD 
and  stored  as  an  M  by  K  array  of  analog  charge  packets. 

By  the  nature  of  charge  transfer  technology,  the  applica¬ 
tion  of  clocking  waveforms  to  the  CCD  electrode  structure 
causes  potential  wells  to  "move"  in  such  a  way  that  these 
charge  packets  are  shifted  along  the  CCD  structure  (upward 
in  the  orientation  of  Figure  I ).  While  the  exact  nature  of 
these  waveforms  varies  from  manufacturer  to  manufacturer 
and  from  one  model  to  the  next,  let  us  assume  that  the  appro¬ 
priate  clocking  pulses  are  provided  so  as  to  shift  these  charge 
packets  vertically  at  a  rate  of  vVCI  CCD  cells/sec.  Let  us  fur¬ 
ther  stipulate  that  vvcr  =  Vjn  .  i.e..  that  the  vertical  shift  rate 
(in  CCD  cells/sec)  is  the  same  as  the  input  signal  sample  rate 
(in  samples/sec).  As  these  charge  packets  travel  upward, 
more  and  more  charge  is  added  to  them  due  to  the  time- 
varying  pattern  illuminating  the  CCD.**  The  net  result  is 
that  each  column  of  the  CCD  performs  a  running  sum  of  the 
products  of  the  time-varying  LED  radiance  and  the  space- 
varying  transmittance  values  of  that  column  of  the  mask.  In 

*Thc  intensity  transmittance  of  the  m.k1*1  cell  of  the  photographic 
transparency  is  proportional  to  the  m.k'h  element  of  the  matrix 
(I'm  g|  This  can  be  achieved  either  by  spatially  modulating  the 
grey  scale  of  the  film  (which  requires  precise  control  of  the  photo¬ 
graphic  exposure  and  development  processes)  or  by  varying  the  area 
of  transparent  apertures  in  an  otherwise  opaque  background  (c.g.. 
spatial  pulse-width  modulation).  For  experimental  case,  the  authors 
chose  the  latter  approach. 

**This  mode  of  operation  of  the  CCD  is  designated  the  "shift  and  add" 
mode  or  “time  delay  and  integration  (TDI)"  mode,  or  even  more  sim¬ 
ply.  the  "processing"  mode  to  differentiate  it  from  the  more  tradi¬ 
tional  imaging  mode. 
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particular,  let  us  consider  the  k  =  nltl  row  of  these  charge 
packets.  As  this  row  travels  from  bottom  (k  =  n=  1 )  to  top 
(k  =  n  =  K),  its  charge  values  take  on  the  form 


^  *n0+n  ^tn.n  ®m.n0 

n=  1 


1,2,3.  — ,M 
--,-2,-1, 0,1, 2,—  . 


(1) 


where  the  unimportant  constant  of  proportionality  has  been 
ignored  for  simplicity.  Performing  the  operation  of  Eq.  ( 1 ) 
in  a  fast,  efficient  manner  is  the  major  accomplishment  of 
the  electro-optical  processor.  We  shall  show  shortly  that 
this  equation  is  the  basis  of  many  mathematical  operations 
frequently  used  in  a  variety  of  signal  processing  applications. 

The  CCD  chip  contains  a  parallel-to-serial  converter 
(physically  just  a  horizontal  linear  CCD  shift  register  with 
M  input  taps)  and  an  on-chip  amplifier  so  that  this  M- 
element  vector  can  be  read  out  of  the  CCD  chip  as  a  time- 
varying  analog  voltage  waveform.  It  is  important  that  the 
shift  rate,  vj,or,  of  this  horizontal  shift  register  be  at  least 
M  times  the  shift  rate  of  the  vertical  CCD  registers  in  order 
that  it  be  empty  by  the  time  the  next  row  of  charge  packets, 
namely  gm,n0+  I  ■  are  ready  to  be  read  out.  Thus,  for  the 
processing  mode  of  CCD  operation. 


vver  ^ 


vhor 

M 


(2) 


This  relationship  between  vvcr  and  V|lor  and  relatedly  the 
need  for  including  the  arbitrary  time  origin  nQ  in  Eq.  ( I ). 
are  of  fundamental  importance  in  understanding  the  opera¬ 
tion  and  capabilities  of  this  processor.  Figure  2  is  a  timing 
diagram  showing  the  relationship  between  the  input  samples 
to  the  LED,  ln0+n  ■  and  the  output  values.  gm  .  from  the 
CCD's  horizontal  shift  register.  Note  that  this  example  de¬ 
picts  the  most  time-efficient  case  in  which  exactly  M  hori¬ 
zontal  shifts  occur  between  every  vertical  shift  (and  conse¬ 
quently  between  every  input  sample).  Thus  the  inequality 
in  Eq.  (2)  is  an  equality.  This  states  that  M  samples  come 
out  of  the  CCD  for  each  input  sample  going  into  the  LED. 


9mnO _ ■■‘I* . ‘I* . lie--- 

1  •••  •••  M  m,n„ 


no=°  no=1  "o=2 


Figure  2.  Timing  relationship  between  input  to  LED  and 
output  from  CCD  for  the  case  of  sliding-window  transfor¬ 
mations.  Overlap  factor  =  (K  -  I  )/K 


With  reference  to  Figure  2,  the  relationship  between 
these  two  time  series  is  as  follows:  Assume  that  the  processor 
is  turned  on  at  some  arbitrary  discrete  time  origin  nQ  =  0  . 
While  the  first  K  input-signal  samples  are  sequentially  modu¬ 
lating  the  LED  and  the  CCD  vertical  registers  are  shifting  up¬ 
ward  K  times,  the  CCD  horizontal  register  is  spewing  out 
(K-  1)M  "garbage”  values.  Immediately  after  the  K,h  input- 
signal  sample,  the  M  values,  gm  q.  outputted  by  the  CCD 
horizontal  register  are  valid  data  -  constituting  the  M  sums- 
of-products  of  the  input  sequence  f|  to  f^  and  the  M  col¬ 
umns  of  the  mask  hm  k  .  Immediately  after  the  K  + Ist  input- 
signal  sample  modulates  the  LED.  the  M  values.  gnl  j  .  out¬ 
putted  by  the  CCD  horizontal  register  again  are  valid  data 
constituting  the  M  sums-of-products  of  the  input  sequence  f-> 
to  fj-  +  |  and  the  M  columns  of  the  mask.  And  so  on.  There 
are  no  more  "garbage”  terms  after  the  initial  start-up  transi¬ 
ent.  Each  sequence  of  M  values  from  the  CCD  horizontal 
register  is  the  sum-of-products  of  the  previous  K  samples  of 
the  input  signal  with  the  M  signals  comprising  the  M  columns 
of  the  mask.  Thus  every  possible  window  of  the  input  signal. 
fn0+n  .  I  <  n  <  K.  is  multiplied  and  summed  with  the  M  ref¬ 
erence  signals  on  the  mask.  In  other  words  the  device  simul¬ 
taneously  generates  the  M  cross-correlation  Jinn  lions  (e.g.. 
the  cross-correlation  value  for  all  values  of  time  delay)  be¬ 
tween  the  input  signal  and  the  M  reference  signals  on  the  mask. 

An  indication  of  the  power  of  this  module  is  that  it  con¬ 
tains  the  equivalent  of  ( I )  an  M  by  K  word  ROM  memory, 
namely  the  matrix  mask.  (2)  a  bank  of  M  by  K  analog  multi¬ 
pliers.  (3)  a  hank  of  M  by  K  analog  adders.  <4|  a  bank  of  \l 
by  K  analog  storage  cells  to  store  the  cumulative  sums,  and 
(5)  a  parallel-to-serial  converter  to  format  these  M  parallel 
data  streams  into  a  single  sequential  output  all  in  the  space 
of  just  a  few  cubic  inches. 

A  restriction  which  the  use  of  incoherent  optics  places  on 
this  processor  is  that,  since  light  intensity  and  intensity  trans¬ 
mittance  can  only  be  positive  (i.e..  non-negative)  quantities, 
the  input  sequence  fn0+n  a,u'  'he  reference  matrix  hm  n 
must  be  positive  quantities.  However,  this  is  not  unduly  re¬ 
strictive  since  bipolar  quantities  can  still  be  processed  by 
either  ( I )  adding  a  bias  to  them  to  make  them  non-negative, 
or  (2)  separating  them  into  their  positive  and  negative  com¬ 
ponents  and  handling  each  of  these  separately.  Similarly 
complex  quantities  must  be  separated  into  their  real  and 
imaginary  parts.  More  detailed  discussions  of  these  tech¬ 
niques  can  be  found  in  References  6  and  7. 

Typical  Applications 
Multi-channel  Cross-Correlation 

As  previously  mentioned  the  electro-optical  processor  can 
be  used  to  compute  the  cross-correlation  functions  between 
the  input  waveform  fn0+n  and  'be  M  signals  hm  n  compris¬ 
ing  the  M  columns  of  the  mask  siimiliuiwimslv  That  is.  it 
compares  all  time  shifts.  nQ.  of  the  incoming  signal  with  a 
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library  of  M  reference  signals  in  parallel.  Such  a  multi¬ 
channel  cross-correlation  operation  is  very  useful,  for  exam¬ 
ple,  in  an  active  radar  system  where  the  M  reference  signals 
might  be  dopplerized  versions  of  the  transmitted  waveform, 
the  input  to  the  LED  might  be  the  received  waveform,  and 
thus  the  locations  m,nQ  of  peaks  in  the  CCD  output  would 
yield  the  range  (proportional  to  nQ)  and  velocity  (propor¬ 
tional  to  m)  of  any  targets.  Also,  such  a  processor  is  useful 
in  a  non-synchronized  communications  receiver  in  which  all 
time  shifts  of  the  incoming  signal  need  to  be  cross-correlated 
with  all  members  of  a  library  of  possible  codewords. 


Multi-channel  FIR  Filtering 


For  any  given  value  of  m,  Eq.  ( 1 )  defines  the  operation  of 

O 

a  shift-invariant  finite  impulse  response  (FIR)  filter.  Thus 
the  electro-optical  processor  can  be  used  as  a  multi-channel 
FIR  filter;  that  is,  each  of  the  M  columns  of  the  mask  can  be 
a  different  FIR  filter  impulse  response.  Each  column  of  the 
processor  could  then  be.  for  example,  a  low-pass,  high-pass, 
or  band-pass  filter  with  each  column  having  a  different  cut¬ 
off  frequency,  bandwidth,  or  center  frequency,  etc. 


General  Linear  Transformation 


For  any  given  value  of  nQ.  Eq.  ( 1 )  also  defines  the  opera¬ 
tion  of  a  general  linear  filter,  i.e.,  a  “black  box"  characterized 
by  an  impulse  response  hm  n  .  This  filter  or  transformer 
need  not  necessarily  be  shift-invariant,  for  which  the  impulse 
response  would  reduce  to  hm_n,  and  is  therefore  able  to 
treat  a  larger  class  of  useful  linear  transformations.  The  fol¬ 
lowing  list  gives  a  few  examples  of  useful  signal  processing 
transformations  which  can  be  performed  by  the  electro- 
optical  processor  together  with  the  appropriate  form  of  the 
required  matrix  mask  hm  : 

Impulse 

Transform  Response 


Convolution 

hm-n 

Cross-correlation 

*’n-m 

Autocorrelation 

f 

n-m 

Cosine  transform 

cos  (2tr  mn) 

Fourier  transform 

exp  (-2tr  imn) 

Laplace  transform 

exp  (-mn) 

Hankel  transform 

2tr  JQ  (2tr  mn)n 

Putting  this  in  another  way,  the  processor  is  a  linear  system 
whose  impulse  response  is  determined  solely  by  the  optical 
transparency  placed  on  the  CCD.  References  4,  5  and  6  show 
some  typical  masks  for  performing  cosine,  Fourier,  and 
Hadamard-Walsh  transforms. 


Reference  to  Eq.  ( 1 )  shows  that  the  electro-optical  proces¬ 
sor  performs  these  transforms  on  the  input  signal  for  every 
consecutive  value  of  nQ.  In  other  words,  it  performs  a  sliding- 
window  transform;  it  performs  the  matrix-vector  multiplica¬ 
tion  operation  of  Eq.  ( I )  on  every  K-point  window  of  the  in¬ 
put  signal.  This  feature  is  critically  important  in  those 


applications  where  the  time  of  arrival  of  the  incoming  signal 
is  unknown,  and  hence  transformations  must  be  performed 
for  all  possible  arrival  times  (or  at  least  with  a  high  degree  of 
window  overlap).  A  modification  to  the  CCD  architecture 
which  would  allow  transformations  with  any  arbitrary  degree 
of  window  overlap,  and  a  consequent  increase  in  input  data 
rate,  will  be  discussed  shortly. 

A  final  comment  on  this  implementation  of  a  linear  trans¬ 
formation  is  that  any  desired  window  function  (e.g..  Ham¬ 
ming.  Poisson,  Gaussian,  Dolph-Tchebyshev,  or  Kaiser- 

Q 

Bessel)  can  be  incorporated  directly  in  the  transparency 
itself  with  no  need  of  any  additional  electronic  or  optical 
components. 

Current  Configuration 

The  most  compact  implementation -of  the  electro-optical 
processor  is  with  the  matrix  mask  in  direct  contact  with  the 
photosensitive  CCD  surface.*  Because  of  the  small  dimen¬ 
sions  of  the  CCD  photosites  (typically  I  mil  by  I  mil)  very 
close  contact  «  I  mil)  is  required  in  order  to  minimize 
smearing  due  to  optical  diffraction.  This,  and  the  necessity 
for  critical  alignment,  strongly  suggest  that  the  mask  should 
be  fabricated  directly  on  the  CCD  surface  itself,  as  a  final 
step  in  the  CCD  fabrication  process.  This  yields  a  CCD  with 
an  integral  pre-aligned  matrix  mask. 

In  1977.  RCA  Laboratories  (David  Samoff  Research 
Center,  Princeton.  New  Jersey)  under  contract  to  the  Naval 
Ocean  Systems  Center  (San  Diego.  California)  fabricated 
masks  onto  RCA's  commercially  available  model  SID  52501 
CCD.  These  devices  have  K  =  5  1 2  by  M  =  320  resolution  cells 
on  1.2  mil  centers.  The  first  step  in  the  process  was  the 
preparation  of  a  photolithographic  mask  containing  the  de¬ 
sired  5  1 2  by  320  element  matrix  of  apertures.  hm  ^ .  The 
remaining  steps  consisted  of  ( I )  depositing  a  protective  oxide 

O 

layer  on  the  CCD  structure.  (2)  depositing  an  opaque  1 200  A 
thick  chromium  film  onto  this  layer,  (3)  imaging  the  photo¬ 
lithographic  mask  onto  this  film  and  etching  away  the  mask 
apertures  using  standard  photoresist  techniques,  and  (4)  eov- 

O 

ering  the  resulting  device  with  a  protective  10.000  A  thick 
SiOs  layer.  The  range  of  aperture  sizes  achieved  in  this  effort 
was  100  discrete  distinguishable  levels  and  improvements  in 
the  technique  are  in  progress  to  push  this  to  about  256  dis¬ 
crete  levels. 

With  such  a  “masked  CCD"  the  only  remaining  compo¬ 
nent  in  Figure  I  is  the  illumination  device.  Successful  results 
have  been  obtained  with  both  the  Monsanto  MV-4  visible 
LED  and  the  Texas  Instruments  TIXL  I  2  P-N  Gallium  Arsen¬ 
ide  infrared-emitting  diode.  The  electro-optical  processor 
geometry  shown  in  Figure  3  has  proven  highly  attractive  since 
(a)  it  maximizes  the  optical  path  from  the  LIT)  to  the  CCD 

•For  those  applications  requiring  programmability  of  the  matrix,  the 
mask  may  take  the  form  of  a  photographic  transparency  imaged  onto 
the  CCD  via  imaging  optics  or  fiber  optics. 
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Figure  3.  Current  configuration  for  the  fixed-mask  electro- 
optical  signal  processor. 

(thereby  maximizing  the  uniformity  of  illumination)  for  a 
given  module  height  by  using  a  folded  optical  system,  (b)  its 
total  volume  is  only  two  cubic  inches,  and  (c)  the  CCD  and 
LED  electrical  contact  pins  are  on  the  same  side,  thereby  per¬ 
mitting  the  unit  to  be  readily  soldered  in  place  on  a  circuit 
hoard. 

This  processor  is  currently  being  tested  and  evaluated. 
Factors  which  adversely  affect  performance  and  which  there¬ 
fore  need  to  be  measured  and  controlled,  are 

LED  linearity. 

spatial  uniformity  of  illumination. 

CCD  spatial  variations  in  sensitivity. 

CCD  charge  transfer  efficiency,  and 
CCD  dynamic  range. 

Expected  performance  parameters  are 

vin  =  Vver  =  32.000  samples/sec 
=  1 0,240.000  shilts/scc 
dynamic  range  >  50  dB 
nonlinearity  <  I'.  . 

Such  a  processing  module  has  a  processing  rate  of  32.000  X 
51  2  X  320  »  5  X  10‘>  multiplications/sec. 

A  Modified  CCD  Architecture 

As  discussed  earlier,  this  electro-optical  processor  multi¬ 
plies  every  possible  time  shift.  n().  of  a  K-point  window  of 
the  input  signal  fn()+n  by  the  M  different  K-point  reference 
signals  hm  ^  on  the  mask  as  described  by  Eq.  ( I ).  This 
"sliding  window"  operation  is  vitally  important  in  the  multi¬ 
channel  cross-correlation  application  but  is  very  wasteful  in 
the  FIR  filtering  and  linear  transformation  applications.  The 
overlap  factor  between  successive  windows  is  given  by 
(K  - 1  )/K  which,  for  K  =  5 1 2.  is  09.8 '3.  Contrast  this  with  the 
5CK/  overlap  factor  typically  used  in  most  Fourier  transform 
operations  where  a  good  window  function  is  used.  Thus,  the 
successive  rows  of  data  gm  nQ  coming  out  of  the  CCD 


contain  highly  redundant  information.  In  other  words,  if  we 
view  a  Fourier  transformer  as  a  bank  of  parallel  band-pass 
filters  we  can  sample  the  filter  outputs  at  a  sample  spacing 
equal  to  the  inverse  of  the  length  of  the  window  function 
used;  whereas  the  electro-optical  processor  samples  them  at 
the  input  signal  sample  rate.  That  is,  the  processor  gives  one 
complete  row  of  Fourier  coefficients  out  for  each  new  sam¬ 
ple  value  read  in. 

It  is  relatively  trivial  to  re-sample  the  coefficients  coming 
out  of  the  CCD,  i.e.,  to  consider  only  every  K/2**1  row  if  a 
507?  overlap  factor  is  desired,  but  this  still  leaves  the  proces¬ 
sor's  input  rate  limited  to  v|)or/M  (i.e.,  to  32.000  samples/sec 
for  the  components  in  the  current  implementation).  The 
bottleneck  in  throughput  rate  is  the  on-chip  parallel-to-serial 
converter.  The  vertical  registers  of  the  CCD  are  capable  of 
shifting  much  faster  than  V|10r/M  shifts/sec  but  are  restricted 
to  this  rate  by  the  time  required  to  read  each  row  of  samples 
out  of  the  horizontal  shift  register.  A  vertical  shift  rate  faster 
than  V|10r/M  would  result  in  charge  packets  from  the  succeed¬ 
ing  row  of  coefficients  being  added  to  the  desired  coefficients 
shifting  across  the  horizontal  register.  Such  contamination 
would  severely  degrade  the  desired  coefficients. 

What  is  needed  is  a  way  of  disearding  the  unwanted  rows 
of  data  coming  out  of  the  tops  of  the  vertical  shift  registers 
without  contaminating  the  desired  row  being  read  out  of  the 
horizontal  register.  This  could  be  achieved  by  incorporating 
a  “controllable  charge  dump"  gate  structure  between  the  ver¬ 
tical  and  horizontal  registers.  This  would  allow  all  unwanted 
output  vectors  g|,1<n  to  be  dumped  to  ground  except  when  a 
predetermined  trigger  is  applied  to  route  the  desired  row  of 
values  into  the  horizontal  register.  To  this  end.  the  Naval 
Ocean  Systems  Center,  under  sponsorship  from  the  Naval 
I  lectronics  Systems  Command.  Code  304.  is  currently  re¬ 
questing  proposals  to  design  and  fabricate  such  a  modified 
CCD  architecture  in  FY  1070, 

For  a  square  CCD  array  (i.e..  M  =  K).  (his  would  allow 
v  (and  consequently  Vjn )  to  be  equal  to  V|)or.  This  case  is 
depicted  in  the  timing  diagram  of  Figure  4  and  corresponds 
to  adjacent-window  transformations  (i.e..  overlap  factor  =  0). 


Figure  4.  Timing  relationship  between  input  to  LED  and 
output  from  CCD  for  the  case  of  adjacent-window  trans¬ 
formations  (i.e.,  M  =  K  and  vver  =  V|)or). 
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This  would  provide  a  compact,  inexpensive  Fourier  transform 
module,  for  example,  capable  of  processing  input  data  rates 
as  high  as  10  MHz.  Any  arbitrary  degree  of  window  overlap 
can  be  provided  for  any  arbitrary  CCD  size  (i.e.,  any  arbitrary 
M  and  K)  by  appropriate  selection  of  vyer  and  vhor  as  shown 
in  the  generalized  timing  diagram  of  Figure  5. 


no=°  nQ  =  A  n0  =  2A 


Figure  5.  Timing  relationship  between  input  to  LED  and 
output  from  CCD  for  the  case  of  transformations  with 
arbitrary  overlapping  of  windows.  Overlap  factor  = 

(K-  M  vver/vhor)/K  =  (vhor-vver)/vhor  for  K  =  M. 

Conclusions 

This  paper  has  qualitatively  described  a  new  signal  proces¬ 
sing  module  capable  of  performing  a  variety  of  useful  opera¬ 
tions  at  high  speeds  and  in  a  compact  package.  However, 
possibly  the  best  feature  of  such  devices  is  their  potential  low 
cost.  As  CCD  yields  increase,  their  prices  will  drop,  thus 
yielding  electro-optical  processors  costing  only  a  few  hundred 
dollars. 
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Introduction 

The  electro-optical  processor  (EOP)  developed 
by  the  authors  is  a  compact,  flexible  device  designed 
for  use  in  many  signal  processing  applications.  Prior 
advances  have  been  described  in  references  1-6.  This 
paper  will  provide  an  overview  of  the  history  and 
mathematical  basis  of  the  device,  present  the  current 
state  of  implementation,  and  indicate  promising  new 
investigation  areas. 

The  EOP  consists  of  an  incoherent  light  source 
(e  g.,  an  LED),  a  photo-generated  mask,  and  an  area- 
atTay  charge  coupled  device  (CCD).  The  package 
which  encloses  these  elements  is  approximately  two 
cubic  inches  in  volume.  Supporting  circuitry  is  con¬ 
tained  on  three  printed  circuit  cards  The  current 
implementation  is  capable  of  performing  5  X  109 
analog  multiplications  per  second. 


Development  History 


During  the  mid-1960’s  the  scientific  literature 
on  optical  signal  processing  emphasised  the  usefulness 
of  coherent  optical  methods  in  performing  various 
mathematical  operations  such  as  two-dimensional 
Fourier  transforms  and  matched  spatial  filtering. 
Investigations  in  this  regard  have  continued,  at  an 
increasing  pace,  with  good  results  achieved.  Some¬ 
what  later,  investigators,  disturbed  by  the  expense, 
size,  and  critical  environmental  controls  usually 
necessary  in  coherent  systems,  recognized  that  in 
some  applications  a  non-coherent  approach  could  be 
substituted.  This  involves  simply  decomposing  all 
complex  bipolar  signals  and  operations  into  their  real 
non-negative  parts,  as  is  done  in  digital  computers, 
thereby  eliminating  the  need  to  operate  directly  with 
complex  quantities. 

Early  in  the  development  of  non-coherent  optical 
processors,  systems  utilizing  television  cameras  and 
line-array  CCD’s  in  conjunction  with  oscillating 
mirrors  were  successfully  demonstrated.*’3,3  The 


availability  of  area-array  CCD’s  operating  in  the  shift- 
and-add  (or  time-delay  and  integration)  mode  permit¬ 
ted  removal  of  the  oscillating  mirrors.  But  even 
with  this  elimination  of  the  mechanical  moving  parts 
these  experimental  systems  still  required  a  significant 
amount  of  space  to  accommodate  the  condensing  and 
imaging  optics  necessary  to  image  the  mask  onto  the 
face  of  the  CCD.  Now,  the  masks  are  fabricated 
directly  onto  the  CCD  obviating  the  need  for  any 
lenses  and  thereby  allowing  non-coherent  processing 
modules  to  be  packaged  in  just  a  few  cubic  inches. 

Mathematical  Basis 

The  generalized  equation  for  a  linear 
transformation 

Jh(m,n)  f(n)  dn  =  g(m)  (1) 

can  be  rewritten  in  its  discrete  form  as 
N 

2  ^m.n  *n  “  ^m 
n=l 

This  can  be  further  stated  as  a  series  of  sums  of 
products, 

N 

2  *M,n  *n  ""  gl  ’ 
n=l 

N 

2  ^2,n  *n  ~  g2  ’ 
n=l 


N 

2  hM,n  fn  =  (3) 

n=l 


•Published  in  the  Proceedings  of  the 
1978  Interna tional  Conference  on  the  Application  of 
Charge  Coupled  Devices.  San  Diego,  C  alifornia. 

25-27  October  1978. 


38 


T 


Some  of  the  forms  of  the  matrix  hm  n  used  to 
perform  various  linear  transformations  are 


Impulse  Response 

^m-n 

^n-m 

*n-m 

cos  (2irmn) 
exp  (-2irimn) 
exp  l-mnl 
2irJ0(2irmn)n 


Transformation 

Convolution 
Cross -correlation 
Autocorrelation 
Cosine  transform 
Fourier  transform 
Laplace  transform 
Hankel  transform 


The  EOP,  as  shown  schematically  in  Figure  1 , 
performs  any  of  these  transformations  as  determined 
by  the  mask  inserted  between  the  LED  and  the  CCD. 
That  is.  each  of  the  M  columns  of  the  mask  and  CCD 
performs  one  of  the  M  sums  of  products  of  Eq.  (3). 
The  following  paragraphs  qualitatively  show  the 
EOP's  operation.  References  4  and  5  contain  a  more 
rigorous  analysis 

In  the  EOP.  the  multiplications  indicated  occur 
as  light  passes  through  a  partially  transparent  mask 
containing  the  matrix  hm  n  as  its  transmittance 
function.  The  intensity  of  the  light  emanating  from 
the  mask  is  indeed  the  product  of  the  irradiance  of 
the  original  light  and  the  intensity  transmittance  of 
the  mask.  Summation  occurs  through  the  ability  of 
the  CCD  to  convert  these  photons  into  charge  packets, 
with  a  direct  linear  relationship  between  total  light 
impinging  on  a  cell  and  the  amount  of  charge  accum¬ 
ulated.  By  shifting  the  charge  packets  vertically 
upward  by  one  resolution  cell  per  input  sample  to  the 


LED,  addition  of  the  proper  terms  through  the 
addition  of  these  charge  packets  is  accomplished. 

To  illustrate  the  principal  of  operation  consider 
the  most  trivial  cases  in  which  M  =  1  (i.e.,  consider 
only  a  single  column  of  the  mask  and  CCD)  and 
N  =  3.  The  light  source  is  consecutively  intensity 
modulated  by  three  input  samples  fj,  fj,  and  f^. 
This  light  illuminates  a  mask  with  transmittance 
values  h| ,  hj,  and  hj  and  falls  upon  the  three  CCD 
cells  C],  C2>  and  c^.  Since  the  incoming  light  floods 
the  mask,  at  the  end  of  the  first  period  the  charge  in 
cells  C|  through  c^ is  proportional  to  f|hj,f|h2,and 
f]hj  respectively.  After  a  shift  of  packets  from  Cj  to 
Co,  C2  to  C3,  and  from  C3  out  to  the  horizontal  shift 
register,  and  another  flood  illumination  by  the 
input  signal  f->,  the  charge  in  cells  C]  through  Cj  is 
proportional  to  fjhj ,  fj  hj  +  f2h->,  and  fjl^  +  f2li3. 
Another  shift  and  illumination  yields  three  charge 
packets  proportional  to  fjhj,  f 2^ |  +  f3h->,  and  fjhj 
+  f2h2  +  f3)>3-  Upon  the  next  shift,  this  last  value  is 
shifted  into  the  horizontal  shift  register  and  read  out 
of  the  CCD.  It  should  be  recognized  as  one  of  the 
sums  of  products  of  Eq.  (3).  By  extension,  it  can  be 
seen  that  an  M  X  N  resolution-cell  CCD,  with  an 
appropriate  mask  can  be  used  to  generate  a  series  of 
M  sums  of  N  products. 

To  continue  the  trivial  example  considered 
above,  if  a  fourth  input  sample  modulates  the  LED 
and  a  fourth  shifting  of  charge  packets  takes  place, 
the  resulting  values  in  cj  through  C3  will  be  propor¬ 
tional  to  f4hj,  f^h |  +  f4h-,,  and  f2hj  +  f^h-)  +  ^3. 
This  last  value,  which  will  be  read*out  of  the  CCD  on 
the  next  vertical  shift,  is  again  one  of  the  sums  of 
products  of  Eq.  (3)  but  with  the  time  variable  of  the 
input  sequence  incremented  by  one.  After  the  initial 
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Figure  1 .  An  Electro-Optical  Processor  (EOP). 
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start-up  transient,  each  new  input  sample  to  the  EOP 
produces  a  set  of  sums  of  products  terms  of  Eq.  (3) 
but  where  the  summation  takes  place  over  the 
previous  N  input  samples  In  other  words,  the  EOP 
performs  sliding-window  transformations.  Or,  in  a 
pattern  recognition  operation,  the  EOP  generates  the 
complete  cross-correlation  function  (i.e.,  the  cross¬ 
correlation  value  for  all  time  delays)  between  the 
input  and  each  member  of  the  library  of  reference 
signatures  on  the  mask. 

It  was  mentioned  earlier  that  signal  decomposi¬ 
tion  was  required  in  those  classes  of  linear  operations 
which  call  for  both  real  and  imaginary  and/or  positive 
and  negative  components  In  the  EOP,  this  decom¬ 
position  is  achieved  by  using  separate  areas  of  the 
mask  for  each  component,  and  by  biasing  the  light 
source,  with  the  results  appropriately  combined  by 
electronic  circuitry  at  the  output  of  the  CCD  6 

Current  State  Of  Development 

In  1977,  RCA  Laboratories  (David  Samoff 
Research  Center,  Princeton,  New  Jersey)  under 
contract  to  the  Naval  Ocean  Systems  Center  (San 


Diego,  California)  fabricated  masks  onto  RCA’s 
commercially  available  model  SID  52501  CCD. 
These  devices  have  N  =  512  by  M  =  320  resolution 
cells  on  1.2  mil  centers.  The  first  step  in  the  process 
was  the  preparation  of  a  photolithographic  mask 
containing  the  desired  512  by  320  element  matrix  of 
apertures  h^  p  The  remaining  steps  consisted  of 
( 1 )  depositing  a  protective  oxide  layer  on  the  CCD 
structure,  (2)  depositing  an  opaque  1200  A  thick 
chromium  film  onto  this  layer,  (3)  imaging  the  photo¬ 
lithographic  mask  onto  this  film  and  etching  away  the 
mask  apertures  using  standard  photoresist  techniques, 
and  (4)  covering  the  resulting  device  with  a  protective 
1 0,000  A  thick  Si02  layer.  The  range  of  aperture  sizes 
achieved  in  this  effort  was  100  discrete  distinguisable 
levels  and  improvements  in  the  technique  are  in 
progress  to  push  this  to  about  256  discrete  levels. 

Figure  2  shows  a  comer  of  the  resulting  “masked 
CCD”.  It  is  interesting  to  note  the  three-phase  elec¬ 
trode  structure  of  the  horizontal  shift  register  at  the 
top  of  the  photograph  and  the  electrode  structure 
(3  electrodes  per  resolution  cell)  of  the  vertical 
registers  showing  through  the  transparent  apertures 
of  the  mask. 
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With  such  a  "masked  CCD"  the  only  remaining 
component  in  Figure  1  is  the  illumination  device. 
Successful  results  have  been  obtained  with  the  Mon¬ 
santo  MV-4  visible  LED.  The  housing  which  holds 
the  LED  and  CCD  also  serves  as  an  optical  cavity  to 
help  provide  a  uniform  illumination  across  the  CCD 
surface  Several  optical  cavity  geometries  are  cur¬ 
rently  being  theoretically  and  experimentally  investi¬ 
gated  with  respect  to  maximizing  uniformity  and 
efficiency  in  a  minimum  size.  As  discussed  more 
fully  in  reference  7,  both  the  cylindrical  and  spherical 
geometries  with  appropriately  placed  baffles  appear 
capable  of  reducing  illumination  nonuniformities 
across  the  CCD  surface  to  under  1%  in  a  housing  size 
of  2  cubic  inches 

This  processor,  with  a  cylindrical  cavity  geom¬ 
etry ,  is  currently  being  tested  and  evaluated.  Factors 
which  adversely  affect  performance  and  which 
therefore  need  to  be  measured  and  controlled,  are 

-  LED  linearity, 

-  spatial  uniformity  of  illumination, 

-  CCD  spatial  variations  in  sensitivity, 

-  optical  crosstalk  between  CCD  cells 

-  CCD  charge  transfer  efficiency,  and 

-  CCD  dynamic  range. 

Expected  performance  parameters  are 

vin  =  vver  =  32.000  samples/sec 

vhor=  '0,240,000  shifts/sec 

dynamic  range  >  50  dB,  and 

nonlinearity  <  1%. 

Such  a  processing  module  has  a  processing  rate  of 
32,000  X  512  X  320  5  X  10^  multiplications/sec. 

The  following  section  discusses  the  external  cir¬ 
cuitry  designed  to  drive  this  EOP  unit.  No  attempt 
has  been  made  to  miniaturize  this  circuitry ;  instead 
the  intention  has  been  to  produce  a  general-purpose 
module  in  which  changes  in  timing  and  operation 


could  be  incorporated  by  simple  rewiring  or  chang¬ 
ing  a  few  components  without  the  need  for  new 
printed  circuit  boards.  This  modular  design  ap¬ 
proach,  in  which  the  circuitry  used  to  operate  the 
EOP  as,  say,  a  Fourier  transformer  is  only  a  slightly 
modified  version  of  that  used  to  operate  the  EOP 
(with  a  different  mask)  as,  say,  a  multi-channel  cross- 
corTelator,  has  saved  considerable  time  in  fabricating 
and  verifying  new  processing  systems. 

An  Overview  of  System  Implementation 

In  the  present  implementation,  the  external 
circuitry  required  to  operate  the  EOP  is  partitioned 
into  three  functional  subsystems.  These  are:  (1) 
master  timing,  (2)  processor  driving,  and  (3)  signal 
conditioning.  As  their  names  imply,  the  master 
timing  circuit  provides  the  required  synchronization 
signals  for  the  overall  system,  the  processor  driver 
controls  the  LED  operation  and  generates  the  nec¬ 
essary  clocking  waveforms  for  the  CCD,  and  the 
signal  conditioner  refines  the  CCD  read-out  to  be 
suitable  for  display  and/or  further  processing.  Each 
subsystem  occupies  one  10  in.  by  4  in.  printed  cir¬ 
cuit  card.  In  the  following  paragraphs,  the  function 
and  design  of  each  card  is  described  in  detail. 

Master  Timing 

There  are  two  functions  for  the  master  timing 
card.  First,  it  defines  the  clocking  intervak  for  CCD 
operation.  Therefore  in  this  respect,  the  design  is 
device  dependent.  The  RCA  CCD  (SIDS 250 1 ) 
employed  requires  both  vertical  and  horizontal  clocks 
to  operate.  The  alternating  intervak  of  vertical  and 
horizontal  clocks  comprise  the  Horizontal  Gate 
as  shown  in  Figure  3.  The  vertical  clock  is  active 
only  between  the  horizontal  clock  intervals. 

Secondly,  the  master  timing  provides  the  mark¬ 
ers  required  to  relate  the  CCD  read-out  to  the  mask 
columns  or  references.  The  CCD  read-out  is  the 
computational  result.  gm-  from  the  EOP:  that  is. 
it  is  a  sequence  of  M-element  vectors  in  the  form  of 
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Figure  3,  Horizontal  Gate. 
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a  time-varying  analog  waveform.  These  markers  en¬ 
able  the  system  to  recognize  which  point  in  time  on 
this  analog  waveform  represents  the  result  of  the 
running  sum  of  products  of  the  input  sequence  fR  and 
the  m'1'  mask  column.  The  numbers  and  kinds  of 
timing  markers  required  are  application  dependent. 
The  card  under  discussion  is  designed  for  general- 
purpose  usage,  in  that  only  minor  wiring  modifications 
are  required  to  provide  the  markers  appropriate  for. 
for  example,  spectrum  analysis,  finite  impulse  response 
I  FIR  i  filtering,  or  cross-correlation  operations. 

This  card  provides  the  following  outputs. 

(I)  Horizontal  Gate,  shown  in  Fig.  3.  This 
defines  the  active  horizontal  clock  interval 
with  a  duration  of  323  horizontal  clock 
pulses  and  the  active  vertical  clock  interval 
with  a  duration  of  7  horizontal  clock  pulses. 
(Each  horizontal  or  vertical  clock  pulse 
results  in  one  shift  along  the  horizontal  or 
vertical  register  respectively.  However,  as 
will  be  discussed  later,  a  vertical  clock  pulse 
occupies  7  horizontal  clock  periods  ) 


black  level  (containing  only  thermally 
generated  “dark  current”)  to  be  used  for 
dc  restoration. 

(5)  Frame  Synchronization  Pulse.  This  occurs 
at  the  start  of  each  5 1 2-line  frame  and  has 
a  duration  of  one  line  (i.e.,  one  period  of 
the  Horizontal  Gate).  A  Frame  Time  is  the 
time  required  to  shift  a  charge  packet  from 
the  bottom  to  the  top  row  of  CCD  cells. 
When  the  electro-optical  processor  is  oper¬ 
ated  in  the  so-called  TV  mode,  for  initial 
check  out  purposes,  the  Frame  Synchroni¬ 
zation  Pulse  is  used  as  the  input  to  the  LED 
so  that  the  LED  flashes  on  only  once 
during  the  entire  frame.  Conveniently,  this 
results  in  an  image  of  the  mask  being  read 
out  from  the  CCD  (i.e.,  each  mask  column 
has  been  correlated  with  a  Dirac  delta  func¬ 
tion  thus  the  correlation  output  is  the  mask 
pattern  itself). 

The  circuitry  to  provide  these  outputs  is  shown 


(2)  Valid  Row  Gate  This  defines  the  location 
of  the  desired  rows  of  the  CCD  output.  In 
many  applications  only  certain  rows  of  out¬ 
puts  are  wanted.  For  example,  in  an  EOP 
performing  a  512-point  Fourier  transform 
with  a  759t  window-overlap  factor,  only 
every  128^'row  of  output  Fourier  coeffi¬ 
cients.  gm.  would  be  desired 

(3)  Valid  Column  Gate.  This  defines  the  loca¬ 
tion  of  the  desired  columns  in  the  CCD 
read-out  line  In  many  applications,  only 

a  specific  sequence  of  columns  are  wanted. 
For  example,  in  a  spectrum  analyzer,  one 
may  be  interested  in  only  a  subset  of  the 
frequency  components,  or  in  a  cross¬ 
correlator.  one  may  wish  to  ignore  the 
columns  used  to  obtain  normalization 
constants  and  zero  levels. 

(4)  Normalization  &  Dark-Level  Pulses.  These 
are  time  signals  indicating  where  the 
normalization  value  and  dark-level  value 
occur  on  a  line  basis.  The  former  is  simply 
the  sum  of  the  previous  N  samples  of  the 
input  signal  and  can  be  obtained  simply  by 
letting  hm  n  =  I  for  one  column  of  the 
mask  This  value  is  necessary  in  normaliz- 
ing  a  cross-correlation  function.  The  dark-j 
level  value  provides  a  convenient  zero  level. 
As  stated  earlier,  there  are  323  horizontal 
clock  pulses,  and  322  states  in  the  hori¬ 
zontal  register  of  the  CCD.  The  one  extra 
clock  pulse  is  for  shifting  out  the  video 


in  Figure  4.  It  consists  of  two  sections,  the  horizon¬ 
tal  clock  section  and  vertical  clock  section  -  each 
composed  of  counters,  decoders  and  combinatorial 
logic. 

The  basic  horizontal  timing  is  maintained  by  a 
set  of  horizontal  counters  They  provide  330  distinct 
clock  periods.  With  this  configuration,  any  of  the 
330  clock  periods  composing  the  line  can  be  pin¬ 
pointed  by  the  use  of  decoders  and  logic  circuits.  In 
this  manner,  the  various  horizontal  timing  signals  are 
generated.  Furthermore,  the  horizontal  section 
produces  a  vertical  pulse  occurring  at  the  line  rate. 
This  pulse  is  used  to  increment  the  line  counters  in 
the  vertical  section  to  maintain  512  distinct  line 
periods.  The  counter  output  is  decoded  to  produce 
the  Frame  Synchronization  Pulse  and  the  Valid  Row 
Gate.  Additionally  a  high  frequency  clock  running 
at  four  times  the  horizontal  clock  frequency  is  avail¬ 
able  for  producing  the  three-phase  horizontal  clocks 
in  the  processor  driver  card. 

Processor  Driver 

The  second  card,  called  the  processor  driver, 
contains:  (1 )  the  LED  modulation  circuitry.  (2)  the 
EOP  module  itself,  and  (3)  the  circuitry  to  generate 
the  CCD  drive  waveforms. 

The  input  signal  to  the  processor  system  is 
applied  to  the  LED  modulator  which  controls  the 
optical  output  of  the  LED.  The  EOP  module  con¬ 
tains  the  LED  and  the  masked  CCD  enclosed  in  an 
optical  cavity,  as  described  in  an  earlier  section.  The 
waveform  circuitry  is  divided  into  the  horizontal  and 
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vertical  clocking  sections  The  former  produces  the 
three-phase,  two-level  horizontal  clocks  and  the  latter 
produces  the  three-phase  multi-level  vertical  clocks.* 

The  circuit  for  the  LED  modulator  and  the 
honzontal  and  vertical  clock  sections  of  the  processor 
driver  are  shown  in  Figure  5.  As  shown  in  the  upper 
right  hand  comer  of  the  circuit  diagram,  the  LED 
driver  consists  of  a  gain  element  (LF357),  a  current 
driver  l  LH0002)  and  the  LED  (Monsanto  MV4  or 
Texas  Instruments  XL12)  connected  in  a  negative 
feedback  fashion  The  LED  current  is  sampled  and 
forced  to  have  the  same  waveshape  as  the  input 
voltage  In  operation,  this  circuit  can  produce  a 
current  pulse  of  0  to  100  ma  through  the  LED  with 
a  rise  time  of  200  nsec. 

The  top  portion  of  Fig.  5  shows  the  horizontal 
clock  section  After  the  horizontal  clocks  are  enabled 
by  the  Honzontal  Gate,  the  three-phase  two-level 
clocking  waveforms  used  to  actually  drive  the  CCD 
are  generated  by  two  stages  of  a  flip-flop  (74LS74) 
connected  in  a  shift-register  fashion.  The  amount  of 
overlap  among  these  phases  is  controlled  by  the 
phasing  of  the  horizontal  clocks  and  the  4X  clocks 
(i.e..  clocks  running  at  four  times  the  horizontal 
clock  rate).  Clock  drivers  (DS0026)  are  used  to 
provide  the  horizontal  drive  capability.  The  outputs 
of  the  clock  drivers  are  capacitively  coupled  to  the 
CCD  so  that  a  dc  bias  to  the  clocks  can  be  inserted. 

The  entire  waveform  can  be  shifted  up  and  down 
with  respect  to  the  CCD  substrate  bias. 

The  lower  portion  of  Fig  5  shows  the  vertical 
clock  section  Each  cycle  of  the  vertical  clocking 
waveform  occupies  the  time  of  seven  horizontal 
clock  periods  Thus  the  vertical  clock  period  is  seven 
times  that  of  the  horizontal  clock  (i.e..  It  takes  seven 
times  longer  to  transfer  all  of  the  charge  packets  in  the 
CCD  vertical  registers  up  by  one  resolution  cell  than  to 
transtcr  those  in  the  horizontal  register  to  the  nght  by 
one  cell  I. 

•The  processor  driver  circuitry  described  herein  is  designed  for  use  with 
the  aforementioned  RCA  SID  52501  CCD.  For  a  more  detailed  discus¬ 
sion  of  the  required  waveforms  than  space  permits  here,  please  refer  to 
the  RCA  specification*  sheet. 


The  Horizontal  Gate  enables  a  binary  counter  to 
maintain  the  seven  district  horizontal  clock  periods 
comprising  each  vertical  clock  cycle.  The  signals 
representing  the  timing  of  each  of  these  distinct 
periods  are  logically  combined  to  control  a  set  of  four 
analog  switches  (RCA  CD4016),  inputs  to  which  are 
fed  with  the  desired  voltage  levels.  One  of  such  wave¬ 
forms  is  illustrated  in  Figure  6. 

For  experimental  purposes,  the  circuitry  to 
generate  the  CCD  drive  waveform  has  been  made  very 
flexible  thereby  allowing  changes  in  various  voltage 
levels  to  be  readily  effected.  Current  amplifiers 
(LH0002)  are  used  as  buffers  between  the  switch  out¬ 
puts  and  the  CCD.  This  circuit  exhibits  rise  and  fall 
times  from  -2  volts  to  +7  volts  of  50  nsec. 

Since  each  card  is  designed  to  be  an  indepen¬ 
dent  unit,  for  modularity  considerations,  a  buffer 
amplifier  (LH0033)  is  used  to  interface  the  CCD  with 
the  signal  conditioning  card. 

Signal  Conditioner 

The  signal  conditioning  card  performs:  (1) 
amplification  of  the  CCD  output,  (2)  dc  restoration. 
(3)  a  sample-and-hold  operation  on  the  amplified 
signal;  and  provides  (4)  output  interfacing  capability. 

The  detail  circuit  and  block  diagrams  are  shown 
in  Figure  7.  A  vido  amplifier  (pA733)  is  used  for 
amplifying  the  CCD  output  from  the  ten  millivolt 
range  to  the  one  volt  range.  Since  the  CCD  output 
is  capacitively  coupled  through  the  amplifier,  the  dc 
level  of  this  resulting  signal  needs  to  be  re-established. 
It  is  convenient  to  utilize  the  dark  level  of  the  CCD 
output  as  a  point  of  reference.  As  mentioned  earlier 
in  the  master  timing  section,  an  extra  horizontal 
clock  is  provided  to  shift  out  the  dark  level.  During 
this  time,  the  dark  level  is  clamped  to  negative  one 
volt  by  use  of  an  analog  switch  (DC,  187)  and  am¬ 
plifier  (LH0033).  The  dc  restored  signal  is  still  con¬ 
taminated  by  feed-through  of  the  clocks,  so  these 
extraneous  signals  must  be  removed  by  using  a  sample- 
and-hold  amplifier  (Date!  SHM-2)  to  obtain  the 
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Figure  6.  One  Phase  of  Vertical  Clock  Waveform. 
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Figure  7.  Signal  Conditioner. 


desired  analog  signal  -  320  sampled  values  per  row. 
Again  a  buffer  amplifier  (LH0033)  is  used  to  inter¬ 
face  the  sample-and-hold  unit  with  any  external 
load  for  output  display  and/or  further  processing. 

New  Investigative  Areas 
A  Modified  CCD  Architecture 

During  the  discussion  of  the  mathematical  basis 
for  the  EOP.  it  was  shown  that  the  EOP  performs 
sliding-window  transformations.  In  those  cases  where 
large  window  overlap  is  desirable,  for  instance  when 
an  incoming  signal  has  some  arbitrary  unknown  start¬ 
ing  point,  this  feature  is  very  useful.  In  those  cases 
where  little  window  overlap  is  desired  some  mechan¬ 
ism  for  flushing  out  the  unwanted  coefficients 
(thereby  increasing  the  bandwidth  capability  of  the 
EOP)  is  desired.  To  this  end,  industry  has  been 
approached  for  solutions,  such  as  a  selectable  gating 
function  inserted  between  the  area  array  and  the 
output  parallel  to  serial  converter.®  If  successful,  the 
vertical  shift  rate  (and  therefore  the  input  sample  rate) 


will  become  much  less  dependent  on  the  horizontal 
shift  rate,  and  in  fact,  limited  only  by  CCD  shift 
drive  capabilities  and  electro-optical  input  signal 
generation. 

A  Real-Time  Programmable  Mask 

Another  promising  investigative  area  is  the  mask. 
It  is  obvious  that  if  a  method  of  generating  masks  in 
real  time  or  near  real  time  existed,  additional  opera¬ 
tions  could  be  performed,  such  as  non-linear  trans¬ 
formations,  recursive  algorithms,  adaptive  solutions, 
and  multi-stimulus  correlations.  Technologies  being 
considered  at  this  time  include  matrix-addressed 
transmissive  liquid  crystal  devices.  Laboratory 
demonstrations  have  been  performed  using  a  Hughes 
liquid  crystal  cell  imaged  onto  the  CCD,  but  con¬ 
siderable  advancement  must  be  made  in  expanding 
the  grey  scale  capability  (i.e .,  dynamic  range),  in 
reducing  the  resolution  cell  size  to  be  compatible 
with  that  of  CCD’s,  and  ultimately  in  fabricating  the 
transmissive  matrix-addressed  mask  directly  on  the 
CCD  surface. 
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Digital  Capability 

The  last  area  to  be  discussed  is  that  of  digital 
operations  with  this  unit  Historically,  optical 
processing  has  been  considered  only  applicable  to 
analog  problems  with  accuracies  corresponding  to 
about  8  bits,  but  work,  is  proceeding  with  mask 
design  and  supporting  algorithms  to  permit  digital 
operations  with  any  desired  degree  of  accuracy. 

Conclusion 

The  EOP  is  a  compact,  low-cost,  low-power, 
high-speed  signal  processing  module  capable  of  per¬ 
forming  a  larger  variety  of  useful  linear  operations 
with  processing  rates  greater  than  10^  multiplications 
per  second  (depending  on  the  array  size  and  maxi¬ 
mum  vertical  shift  rate  of  the  CCD  used).  General- 


purpose  modules,  each  consisting  of  a  two  cubic 
inch  EOP  and  three  cards  of  external  circuitry,  have 
been  fabricated  and  are  now  undergoing  test  and 
evaluation  in  a  variety  of  application  areas  Also, 
further  development  work  is  planned  to  increase 
the  EOP’s  input  rate  from  the  present  32  kHz  to  the 
MHz  range,  to  provide  real-time  programmability  of 
the  mask,  and  to  allow  the  performance  of  digital 
operations  with  arbitrary  accuracy. 
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QUAD  CHIRPED  TRANSVERSAL  FILTER/R5601 


KEY  FEATURES 

Provides  the  real  and  imaginary  convolutions  necessary  for  the  chirp  Z 
transform  algorithm. 

Two  512-stage  charge-coupled  devices. 

Four  mask-programmed  transversal  filters. 

Balanced  differential  outputs. 

Dynamic  range  (peak  signal  to  rms  noise)  of  60  db. 

Sampling  rates  4  KHz  to  2  MHz. 

Filter  weighting  coefficient  accuracy  of  eight  bits  plus  sign. 

22-pin  dual-in-line  package. 


INTRODUCTION 


The  R5601  is  an  MOS  integrated  circuit  which  can  be  used  to  perform  the  bulk  of  the 
computation  required  in  the  calculation  of  a  512-point  Discrete  Fourier  Transform  (DFT).* 

This  circuit  contains  two  separate  512-stage  charge-coupled  devices  which  are  used  to 
implement  four  transversal  filters  using  the  split-electrode  weighting  technique.2 
The  filter  weighting  coefficients  and  internal  circuit  connections  are  configured  so  that  this 
device,  with  additional  off-chip  components,  can  implement  the  chirp  Z  transform  (CZT) 
algorithm  to  calculate  a  512  point  DFT.2 

There  are  two  versions  of  the  R5601  which  are  available.  They  differ  in  the  type  of  window 
which  is  effectively  applied  to  the  data  to  be  transformed.  The  R5601-1  uses  a  rectangular 
window  for  maximum  resolution  in  the  transform  (or  frequency)  domain  with  the  disadvantage 
of  relatively  high  sidelobe  response  (frequency  leaks).3  The  R5601-2  uses  a  Hanning 
window  in  order  to  decrease  these  sidelobes.  However,  this  window  also  decreases  the 
frequency  resolution  of  the  resulting  DFT  by  about  a  factor  of  two.  In  general,  the  R5601-2 
is  intended  for  use  in  spectral  analysis  applications  and  the  R5601-1  for  direct  calculation 
of  the  complex  Fourier  coefficients. 

The  use  of  analog  discrete  time  techniques  to  calculate  Fourier  spectra  allows  an  enormous 
reduction  in  hardware  complexity  and,  therefore,  cost  when  compared  to  all-digital  techniques. 
This  approach  should  find  application  whenever  low  cost,  reduced  size,  high  speed,  or  low 
power  requirements  arc  important. 
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QUAD  CHIRPED  TRANSVERSAL  FILTER/R5601 


DEVICE  DESCRIPTION 


The  R5(i01  is  available  in  a  22  -pin  DIP  package  with  connections  as  shown  in  Figure  1. 

The  functions  of  the  pins  are  given  in  Table  I  and  a  simplified  block  diagram  in  Figure  2. 

The  timing  diagram  is  shown  in  Figure  3  and  a  typical  connection  arrangement  in  Figure  4. 
Figures  5  and  (>  depict  the  device  geometry  and  the  I/O  circuits.  The  chip  contains  two 
CCD  channels  l>ing  adjacent  to  each  other.  Figure  5  is  an  edge  view  of  the  chip  which  has 
been  cut  along  the  length  of  the  channels  to  exhibit  the  four-phase  clock  gates  which  are 
driven  by  an  external  clock  to  transfer  the  signal  charge.  The  input  signal  gates,  ISGR  and 
I  SGI  »  for  real  and  imaginary  channels,  respectively,  are  clocked  with  the  same  phase 
as  02:  therefore,  the  signal  is  sampled  into  both  channels  when  0(9  goes  low.  These  discrete 
signal  packets  then  pass  under  the  input  receiving  gates,  IRG^  and  IRqj,  which  are  the 
receiving  gates  for  the  real  and  imaginary  channels  respectively,  then  sequentially  transfer 
down  the  channel  through  the  potential  wells  which  are  formed  under  each  of  the  four  clock 
phases,  /p  0 9,  0^  and  0^,  through  clocks  applied  from  an  external  source.  After  trans¬ 
ferring  through  the  channel,  the  signal  terminates  into  V  .f|  through  the  output  gate,  VoG 
and  the  switch  labeled  RESET.  The  delay  line  outputs,  labeled  VOI  and  VOR,  are  required 
only  during  testing  and  initial  adjustment. 

Figure  6  depicts  the  top  view  of  a  single  CCD  channel  and  demonstrates  the  split-electrode 
technique  which  implements  the  two  different  filters  on  each  channel.  The  weighting  co¬ 
efficients  of  these  filters  are  determined  by  each  major  gate  electrode  which  is  split  with 
weighting  proportional  to  '.he  difference  in  the  lengths.  The  chirp  waveform  in  dotted  lines 
emphasizes  the  two  separate  weighting  functions  on  the  same  channel,  one  weighting  effective 
during  tf^'s  high  state,  and  the  other  effective  during  0^'s  high  state. 

In  order  to  perform  the  Fourier  transform,  it  is  necessary  to  multiply  the  input  signal  by  a 
complex  chirp  waveform,  convolve  the  real  and  imaginary  parts  with  a  second  complex  chirp 
waveform  and  then  post  multiply  the  output  by  the  same  chirp  waveform  used  in  the  pre¬ 
multiplication.  The  R5601  performs  the  convolution  of  the  real  and  imaginary  input  signals 
with  a  complex  chirp  e^1 n  Vn,  1  -  n  -1  512  for  the  R5601-1  .  .  If  only  the  magnitude  of 
the  spectrum  is  required,  then  the  filter  outputs  are  squared  and  summed.  This  is  true  for 
both  the  R5601-1  and  R5601-2.  Block  diagrams  for  calculating  the  power  spectra  |  Fjq  2  and 
the  complex  Fourier  coefficients  F^  are  shown  in  Figure  7.  The  major  computational 
task  —  the  complex  convolution  —  is  performed  by  the  R5601  convolution  filter  shown  within 
the  dash-line  rectangle. 

To  perform  the  convolution,  the  R5(i01  device  contains  four  512-stage  CCD  mask-programmed 
transversal  filters  which  are  implemented  using  two  CCD  delay  lines.  The  two  different 
filter  functions  on  each  CCD  are  programmed  by  splitting  the  0^  and  0%  clock  electrodes  which 
are  labeled  (0^+ ,  0^~  and  0^' ,  0,~)  as  shown  in  Figure  0.  The  effective  multiplying  factor, 
called  the  weighting  coefficient,  is  proportional  to  the  difference  in  size  of  the  two  portions 
of  a  particular  electrode.  These  weighting  coefficients  are  programmed  to  an  accuracy  of 
eight  bits  plus  sign  (i.e. ,  512  possible  levels).  Since  the  value  of  the  weights  are  determined 
by  the  pattern  of  one  mask  used  in  the  semi-conductor  fabrication  process,  it  is  expected  that 
this  accuracy  level  will  be  independent  of  minor  fabrication  variations  and  will  be  easily 
reproducible. 
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The  sampled  analog  signal  moves  down  the  analog  delay  line  under  control  of  the  four 
clock  phases,  0.,  02,  03  and  04  (See  Figure  6).  The  04  clock  drive  should  be  capacitively 
coupled  by  two  capacitors  to  the  programmed  split  electrodes  (Pin  G),  and  04  (Pin  4) 
which  are  also  the  output  signal  lines.  This  is  also  true  for  the  03  clock  and  the  split 
electrodes  0„+  and  03”.  The  values  of  the  external  capacitors  are  500  pf.  Pins  8,  11,  and 
15  also  require  the  03  clock  waveform,  but  it  should  come  from  a  different  driver  device 
than  the  one  used  for  03+  to  insure  isolation  of  the  output  from  the  input. 

The  outputs  from  the  split  electrodes  (i.e.,  0j+  and  04  ,  03  and  0^  )  are  appropriately 
summed  together  on  the  chip.  The  output  from  similar  clock  electrodes,  for  example 
0  +,  of  the  two  different  delay  lines  are  brought  to  the  same  point  on  the  chip  (See  Figure  6). 
This  connection  performs  the  summing  of  the  signals.  In  the  typical  configuration,  the 
output  from  0  +  and  01“  are  differentially  combined  and  then  sampled  and  held  until  03  data 
is  valid.  The103H  and  03”  are  also  differentially  combined  to  produce  the  other  output  signal. 
Two  more  clock  signals,  02  and  04,  are  minor  phases  required  to  drive  these  devices 
(See  Figure  6).  Two  off-chip  differential  amplifiers  with  a  gain  of  approximately  10  for 
the  5601-1  or  20  for  the  5601-2  are  required  to  obtain  the  output  signals  of  five  volts 
maximum  amplitude. 

The  analog  input  signals  into  the  CCD  device  are  superposed  on  an  approximately  6-volt 
dc  bias  level  and  can  swing  approximately  four  volts  peak-to-peak.  The  input  circuitry 
is  designed  to  be  compatible  with  the  split-electrode  output  techniques  to  give  maximum 
linearity  and  dynamic  range  for  this  device.  In  particular,  the  input  structure  is  designed 
to  compensate  for  the  non-linear  effects  of  the  depletion  capacitance  under  the  CCD  gates 
of  the  filter. 

Pins  10  and  13  are  test  points.  Each  should  be  connected  to  a  lOKohm  resistor  which  is 
connected  to  the  substrate  voltage.  These  pins  provide  the  delayed  value  of  the  input  signal 
during  testing  procedures. 

THE  TRANSFORM  ALGORITHM 

The  complete  Discrete  Fourier  Transform  (DFT)  is  given  by  the  formula 

Fk  =  ''I' nk/N  k  =  0,1,2 _ N-l  (1) 

n-0 

where  either  or  both  Fk  and  fn  may  be  complex.  The  factor  2nk  in  the  exponent  can  be 
replaced  by  its  seemingly  more  complicated  equivalents 


This  substitution  changes  Eq.  1  to  Eq.  3  below,  where  the  pertinent  factors  have  been 
segregated  to  emphasize  the  important  operations. 

_  -ilT  k2/N  N-l  ,,  -iiT  n2/N,  iir(k-n)2/N 
Fk  =  e  (f n  e  >  e 


'there  are  three  operations  indicated  by  Eq.  3 

1.  Multiply  each  corresponding  term  of  the  input  discrete-time  series,  fn, 
by  the  complex  factor,  exponential  (-iTTn2/N);  that  is,  by  the  pair  of 
factors  cos  frn2/N  and  i  sin  rtn^/N,  to  produce  the  combination  term 
in  parentheses.  Let  this  new  complex  sequence  be  called  gn> 

2.  Perf^rjn  a  discrete  convolution  between  the  sequence  gn  and  the  sequence 
eiTrn  ;  this  is  the  portion  within  the  summation. 


3.  Multiply  the  resulting  output  sequence  by  the  final  factor  e 
for  each  point  of  F^. 


-i  tt  k2/N 


These  three  operations  are  indicated  in  the  block  diagram  of  Figure  7,  which  shows 
the  convolution  portion  within  the  dash-line  rectangle.  For  spectral  density,  outputs 
from  the  convolution  filter  are  squared  and  summed;  for  obtaining  the  complete  Fourier 
coefficients,  the  post  multiplier  shown  at  the  right  in  Figure  7  must  be  provided  in  place 
of  the  squaring  function.  For  a  512-point  transform,  the  convolution  portion  requires 
a  minimum  of  2048  multiplications  and  2048  additions  for  each  Fk;  it  thus  represents  the 
bulk  of  the  computational  task.  It  is  this  portion  which  is  performed  by  the  R5601  convolution 
filter. 

Note  that  complex  multiplication  and  complex  convolution  are  linear  operations;  that  is, 
operation  on  a  sum  is  the  same  as  the  sum  of  the  results  of  separate  operations  on  the 
members  of  the  sum.  Thus,  these  operations  can  be  broken  into  corresponding  real  and 
imaginary  (quadrature)  components.  Note  further  that  these  operations  can  be  handled  by 
real  operations,  but  in  separate  real  and  quadrature  (imaginary)  channels.  Thus 

x  .  y  =  (xR  +  ixj)  .  (yR  +  iy^  =  (xRyR  -  xiyi)  +  i  (xRyT  +  xTyR) 

x  *  y  =  (xR  *  yR  -  xr  *  yj)  +  i  (xR  *  yj  +  Xj  *  yR)  and  (3a) 

(x  .  y)  *  z  =  [(xRyR  -  Xjyj)  *  zR  -  (xRyj  +  XjyR)  *  zj  +  l  [(xRyR  -  xft)  *  Zj 

+  <xRyi +  w  *  zr]  (3b) 

These  are  the  operations  indicated  in  Figure  7  and  Eq.  3  where  x  is  the  signal,  y  is  the 
chirp  multiplier  and  z  is  the  convolution  portion  supplied  by  the  R5601. 


A  general,  complex  input,  which  has  separate  real  and  imaginary  (quadrature)  inputs, 
uses  the  four  premultipliers  shown  at  the  left  in  Figure  6,  and  a  single  (complex)  Fourier 


coefficient  is  obtained  for  each  frequency  component  of  the  input.  The  usable  input  frequency 
band  is  from  zero  to  the  sample  frequency,  without  ambiguity,  because  each  input  sample 
has  two  components  (that  is,  there  are  effectively  2  f  sample  values  per  second,  and  the 
Nyquist  criterion  thus  allows  a  band  to  fg  without  aliasing).  Another  way  of  viewing  these 
results  is  to  observe  that  the  input  premultiplier,  with  real  and  quadrature  inputs,  forms  a 
single- sideband  chirp  modulator,  so  that  only  a  single  frequency  occurs  in  the  chirp 
input  to  the  filter  for  each  component  of  the  signal. 


In  most  cases,  however,  inputs  are  real,  the  imaginary  input  is  identically  zero,  and  two 
of  the  multipliers  may  be  deleted.  This  (real)  signal  is  simply  split  into  two  paths.  The 
signal  in  one  path  is  multiplied  by  a  cosine  chirp  (cosir  n2/N)  and  fed  to  the  real  channel  of 
the  filter  input;  the  signal  in  the  second  path  is  multiplied  by  a  negative  sine  chirp  (-sin-rrn2/N) 
and  fed- to  the  imaginary  (quadrature)  channel.  Each  real  input  component,  however,  is 
composed  of  two  complex  components,  one  with  positive  frequency  and  one  negative 
frequency.  For  example,  a  real  signal,  fn  =  cos  2rrnT,  is  composed  of  complex 
components,  1/2  r  and  1/2  J-rrr|l'j  wh0se  imaginary  portions  cancel.  These  two 
components  generate  sum  and  difference  chirp  frequencies,  and  for  every  real  signal  in  the 
band  0  to  f  sample  /  2  the  negative-frequency  component  turns  up  in  the  bano  isample  to 
^sample  /  2  as  an  alias,  thus  restricting  the  useful  input  band  to  components  lying  below 
this  Nyquist  limit  of  fsample  /2. 


The  operations  above  will  perform  the  complete  DFT,  but  to  do  so  requires  careful 
attention  to  the  total  operation.  In  many  cases,  however,  only  the  magnitude  of 
the  power  spectrum  is  required  (no  phase  information),  and  considerable  simplification  is 
possible.  First,  only  the  squared  magnitude  of  the  various  F^  is  required,  so 


2 


(fn 


^-irr  n2/Nj  eirr  (k-n)  2/N 


(4) 


2 

The  final  phase  multiplier  term,  e  k  ^N,may  be  deleted  because  it  has  unit  magnitude 
and  so  does  not  affect  the  amplitude.  Second,  it  is  possible  to  operate  in  a  continuous  fashion 
with  only  slight  approximation,  because  so  many  points  are  involved.  We  step  the  input 
data  each  time  a  new  spectral  component  is  calculated.  For  a  periodic  waveform,  the 
principal  effect  is  a  slight  phase  factor  which  has  only  slight  effect  on  the  result.  Equation  (4) 
then  becomes: 


I  N- 1 


rS 

Fk 


-  £ 


(f. 


n+k 


rr  n2/N.  irf(k-n)2/N 
)  e 


n=0 


(5) 


where  the  designation  FkS  indicates  a  sliding  (or  continuous)  chirp  Z  transform  (2\  The 
chirp-Z  designation  arises  because  of  the  frequency  factors  e* TT*  n  /N  which  effectively 
change  frequency  linearly  through  the  interval  Ot  n<N.  The  phase  and  frequency  are, 
respectively:  ^ 

I  9  =  TTn  /N  radians 

j^n|  =  j^/dn  |=  2Tf  n/N  radians/step 
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For  real  inputs  then,  only  one  (common)  input  is  used.  Further,  each  filter  pair  gives 
one  component  of  a  complex  number,  so  for  the  spectrum  analyzer  implementation,  each 
component  is  squared,  then  added  to  the  other  to  give  the  overall  squared  magnitude.  The 
interior  complex  operations  of  the  convolution  are  separated  into  their  real  and  imaginary 
components  and  handled  separately  in  the  real  and  imaginary  filter  channels.  This  is  the 
situation  addressed  by  the  RC5G01  evaluation  module  which  is  available  from  Reticon. 

The  common  input  f(n)  is  simply  split  into  two  paths  for  processing.  The  signal  in  one 
path  is  multiplied  by  a  cosine  chirp  and  delivered  to  one  pair  of  ,convolver-filters  to 
constitute  the  "real"  channel;  the  signal  in  the  other  path  is  multiplied  by  a  sine  chirp  and 
delivered  to  the  second  pair  of  convolver-filters  to  constitute  the  ^'imaginary"  channel. 

The  outputs  are  combined  on  an  approximate  rms  basis  to  give  the  spectral  density. 

To  see  how  the  convolver-filters  work,  consider  the  chirped  waveform  of  Figure  8  where 
the  cosine  chirp  is  plotted  around  the  circumference  of  a  circle  to  emphasize  the  repetitive 
nature  of  the  chirp. 

Only  the  cosine  chirp  is  shown  in  the  plot  for  clarity.  Its  equation  is 

r  =  2.  5  +  0.  5  cos  Tf  n2/36 
cos 

where  n  =  0,  1,  2. . .  .36  and  where  r  =  2.  5  is  the  "axis"  of  the  waveform.  (For  the 
figure,  N  =  36,  to  keep  the  figure  simple;  in  the  R5601,  N  =  512).  Note  that  the  waveform 
is  accurate  for  0i  n  <  N/2,  and  thereafter  aliasing  makes  the  discrete  points  also  lie 
on  a  mirror  image  of  the  first  18  points,  as  shown  dotted,  whereas  the  continuous  chirp 
curve  has  constantly  increasing  frequency  with  many  more  alternations,  as  shown  by  the 
solid  waveform. 

Now  consider  a  dc  input.  Multiplication  by  the  input  chirp  weighting  produces  time  samples 
exactly  like  those  of  Figure  8.  These  samples  are  fed  into  the  delay  line  of  the  convolver- 
filters.  At  the  end  of  512  clock  periods,  the  samples  in  the  filter  exactly  match  the  filter, 
and  the  output  is  maximum  at  this  instant  in  time  (the  zero-freqie ncy  position). 

Next,  consider  a  fixed  low-frequency  cosine  input  whose  samples  are  f(t)  =  cos  Tf  mn/N, 
with  m  an  arbitrary  (small)  factor.  This  input  is  multiplied,  point  by  point,  with  the  cosine 
and  sine  chirps  to  give  corresponding  outputs  to  the  filter  channels  of 

P  =  cosfrr  mn/N )cos(rf  n2/N) 

=  1/2  cos  ^  (n2  +  mn)  +  cos  (n2  -  mn) 

Q  =  cos(rr  mn/N)  sin^r  n^l) 

=  1/2  sin  — -(n2  +  mn)  +  sin  ^  (n2  -  mn) 

selves  are  chirped  sum  and  difference  frequencies  just  like  those  in 
•d  forward  or  backward  by  the  shift  in  frequency.  Thus,  we  see 


that  there  will  be  a  new  time  slot,  relatively  earlier  or  later  than  that  for  dc,  where  the 
samples  exactly  match  the  filter  to  give  maximum  output.  At  any 

other  position  more  than  a  small  fraction  of  a  time  slot  away,  the  match  grows  rapidly 
worse  to  give  essentially  zero  output  when  all  the  weighted  samples  are  added.  Note  that 
there  is  a  pair  of  responses  from  the  filter,  corresponding  to  the  sum  and  difference  fre¬ 
quencies,  since  forward  rotation  corresponds  to  the  sum,  and  rearward  rotation  corresponds 
to  the  difference,  and  matches  occur  due  to  the  symmetry  of  the  (sampled)  pattern.  As 
the  input  frequency  is  increased,  the  two  responses  move  toward  each  other  until  they  meet 
for  an  input  at  the  Nyquist  frequency  of  fsamPie/2-  There  are  thus  256  time  slots  where 
discrete  frequency  components  can  have  convolution  maxima,  with  frequency  components 
distributed  between  dc  and  the  Nyquist  frequency.  The  resolution  thus  is  dependent  on  the 
number  of  points  in  the  matched  filters.  So  if  we  have  512  points  and  100  KHz  sample 
frequency,  we  would  end  up  with  approximately  200Hz  resolution;  with  10  KHz  we  would 
have  20  cycles  per  sample  time  resolution.  However,  it  should  be  noted  that  a  frequency 
midway  between  adjacent  resolution  cells  spreads  its  energy  equally  between  the  cells; 
in  either  event  the  area  under  the  response  tends  to  remain  constant.  Windowing,  and 
inputs  nonsynchro  nous  with  the  sample  rate  also  tend  to  spread  the  spectrum  and  cause  it 
to  vary  in  amplitude  with  time. 


Why  are  there  four  matched  filters?  Equation  (3b)  indicates  why.  Even  if  the  input  is 
purely  real  (i.e.,  Xj  =  0),  the  correlation  process  still  requires  all  four  components  of 
z;  two  in  the  real  channel  and  two  in  the  imaginary  channel.  Qualitatively,  for  example, 
while  a  real,  even  input  correlates  primarily  with  the  cosine  chirp,  the  sidelobes  do  not 
cancel  out  nearly  as  well  without  the  other  filters  present.  For  more  general  signals, 
each  filter  section  contributes  to  the  main  correlation  as  well. 


SPECIFICATIONS  AND  PERFORMANCE 

The  diagram  of  Figure  4  indicates,  in  part,  the  configuration  and  drive  requirements  of  the 
R560I.  Typical  operating  conditions  are  indicated  in  Table  II.  Input  impedances  to  the 
various  gates  are  principally  capacitive  in  nature.  The  delay-line  outputs  are  connections 
to  active  FET  sources  of  the  output  source  followers,  so  these  outputs  normally  are  connected 
through  load  resistances  (<^10K  ohm)  to  the  substrate  potential.  These  delay-line  outputs 
are  normally  used  only  during  testing  and  initial  adjustment.  The  delayed  signal  appears 
superposed  on  clock  pedestals  as  in  Figure  9  ,  where  the  cross-hatched  portion  represents 
variable  values,  depending  on  the  particular  signal  sample.  If  the  viewing  oscilloscope 
sweep  is  synchronized  to  an  input  sinusoidal  signal,  the  output  should  show  a  stair-step 
approximation  to  the  input  when  biases,  etc.,  are  correctly  adjusted. 


EVALUATION  MODULE 

A  printed-circuit  evaluation  module  is  available  from  Reticon.  This  module  provides  the 
external  functions  of  clock  drive,  premultipliers,  output  buffers,  squaring  circuits,  etc. , 
to  implement  a  spectral-density  evaluator  for  real  inputs.  It  is  self-contained  except  for 
power  supplies.  Its  design  is  a  compromise  between  performance  and  cost,  making  it 
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useful  principally  for  initial  evaluation  and  simple  spectral  analysis.  Its  nominal  sample 
rate  is  100  KHz  but  an  external  input  permits  sampling  at  lower  rates  if  desired.  Further 
details  may  be  found  In  specification  data  for  the  RC-5601  Evaluation  Board. 
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Table  I 

H5601  Pin  Functions 


Pin  No. 

Designation 

Function 

1 

Vsub 

Substrate 

2 

04 

Minor  clock  for  CCD  transport 

3 

COM 

Common  (ground)  connection 

4 

*r 

(-)  weights  of  0j  split-electrode  gates 

5 

0rl 

Reset  gate  for  0^  split  electrodes 

6 

*1+ 

(+)  weights  of  split-electrode  gates 

7 

^GI 

Input  sampling  gate,  imaginary  channel 

8 

trgi 

Input  receiving  gate,  imaginary  channel 

9 

Input  Imaginary 

Signal  input,  imaginary  channel 

10 

VOI 

Delay-line  output,  imaginary  channel 

11 

Reset 

Delay-line  output  reset  control 

12 

VOG 

Delay-line-output- control- gate  bias 

13 

VOR 

Delay-line  output,  real  channel 

14 

Input  Real 

Signal  input,  real  channel 

15 

JRGR 

Input  receiving  gate,  real  channel 

16 

ISGR 

Input  sampling  gate,  real  channel 

17 

*3+ 

(+)  weights  of  03  split-electrode  gates 

18 

0r3 

Reset  gate  for  03  split  electrodes 

19 

<*3~ 

(-)  weights  of  03  split-electrode  gates 

20 

Vdd 

Drain  supply  for  source-follower  buffers 

21 

02 

Minor  clock  for  CCD  transport 

22 

VCB 

Delay-line  corner-control-gate  bias 
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Table  II 

Definitions  and  Electrical  Specification 


Pin 

No. 

Functions 

Symbol 

Parameters 

Min  TVP.  Max. 

Units 

Note 

1 

Substrate 

Sub 

0 

-1.2 

-7 

volts  dc 

3,5 

2 

Minor  phase  clock 

*4 

5 

6 

7 

p-p  volts 

1.4 

3 

Common 

Com 

— 

— 

— 

— 

2 

4 

Major  phase  clock  & 

Signal  Output  Line  for 

K 

10 

12 

14 

p-p  volts 

1.4 

Minus  Signal 

5 

Reset  gate  for  jL  clock 

line 

^rl 

5 

15 

22 

p-p  volts 

1.4 

6 

Major  phase  clock  and 

Signal  Output  Line  for 

K 

10 

12 

14 

p-p  volts 

1,4 

Fbsitive  Signal 

7 

Input  sampling  gate. 

Imaginary 

tsgi 

10 

15 

22 

p-p  volts 

1.4 

8 

Input  receiving  gate. 

Imaginary 

trgi 

0,5ISGI 

0,75ISGI  rSGI 

volts  d.c. 

5 

9 

Signal  input, Imaginary 

IN-Im 

2 

4 

7 

volts  d.c.  (bias) 

10 

Delayed  output.  Imaginary 

VOI 

see  text 

11 

Reset  input  to  Channel 

RESET 

10 

15 

22 

p-p  volts 

1 

Output  Amplifier 

12 

Signal  Channel  Output 

Gate 

VOG 

2 

3 

5 

volts  d.c. 

5 

13 

Delayed  Output,  Real 

VOR 

see  text 

14 

Signal  Input,  Real 

IN-Re 

2 

4 

7 

volts  d.c.  (bias) 

15 

Receiving  Gate,  Real 

1  RGR 

0i5ISGR 

°'75ISGRISGR 

volts  d.c. 

5 

16 

Input  Sampling  Gate, 

Real 

!sgr 

10 

15 

22 

p-p  volts 

1 

17 

Major  Phase  Clock  & 

Signal  Output  Line  for 

*3+ 

10 

12 

14 

p-p  volts 

1.4 

Positive  Signal 

18 

Reset  clock  for  0  clock 

line 

^r3 

5 

15 

22 

p-p  volts 

1,4 

1 
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Table  II  (continued) 

Definitions  and  Electrical  Specification 


pin  Parameters 

No.  Functions _ Symbol  Min  Typ  Max _ Units _ Note 


19 

Major  Phase  Clock  &  Signal  ^  - 
Output  Line  for  Minus  Sigal  3 

10 

12 

14 

p-p  volts 

1,4 

20 

Supply  voltage  for  output  Vdd 

8 

15 

22 

volt  d.  c. 

1,5 

21 

Minor  phase  clock  #2 

5 

6 

7 

p-p  volts 

1 

22 

Comer  bias  VCg 

2 

3 

5 

volts  d.c. 

1,5 

Note:  1.  All  clocks  swing  from  0  to  0.4  volts  low  to  a  high  as  specified. 

2.  Common  is  reference  level  for  all  voltage  measurements  and  is  normally  grounded. 


3.  The  substrate  must  be  at  lowest  potential,  normally  -1.2  volts. 

4.  See  figure  for  clock  waveform  and  timing;  note  that  all  waveforms  must  be 
free  from  over  and  under  shoots  with  edges  as  clean  as  possible. 


5.  These  are  bias  input  nodes  and  must  be  well  by-passed. 
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Table  II  (cont. ) 

Specifications  (25°C) 

INPUTS 

Signal  Inputs 

Symbol 

Typical 

Units 

Bias  Level 

Input  Real 

6 

volts 

Signal  Level 

Input  Imaginary 

4 

volts  (p-p) 

Input  Capacitance 

5 

pf 

Clock  &  Drive 

Split- Electrode  Clock  Amplitude 

30 

volts 

(Applied  to  Capacitors,  see  Fig.  2 

Electrode  Clock  Amplitude 

#2» 

3-5 

volts 

Reset  Clock  Amplitude 

*Rl'*R3Reset 

12-15 

volts 

Clock  Line  Capacitance 

02» 

200 

pf 

+  ± 

Clock/Sense  Line  Capacitance 

I— ‘ 

CO 

500 

pf 

Minimum  Sample  Rate 

4 

KHz 

Maximum  Sample  Rate 

2 

MHz 

Drain  Supply 

vdd 

15 

volts 

Output  Gate 

V0G 

3 

volts 

Comer  Gate 

VCG 

3 

volts 

Substrate  Bias 

^sub 

-1.2 

volts 

Power  Dissipation 

— 

.5 

watts 

OUTPUTS  f01+,  0,~,  0O~, 

1  Signal  Level  (Superimposed  on  Clock  Level) 

.5 

volts  (p-p) 

PERFORMANCE 

• 

• 

Dynamic  Range  (Peak  signal  to  RMS  noise) 

60 

db 

Linearity  (Total  harmonic  distortion 

1 

% 

Tap  weight  Accuracy 

8  bits  plus  sign 

t 

The  Weighting  coefficients  for  the  two  devices  are: 

R5601-1  COS.  C  7T  (n-1)2  ]  Tap  Number  =  n 

SIN  |  512  J  1  ^  n  ±  512 

R5601-2  COS.  f  7T(n-257)2  }  f  1/2  |~1  -  COS  /  2-rT  (n-l)\1 1  1  ^  n  *  512 

SIN  512  j  [  L  V  511  {JJ 


. 
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COMMON  |_  3 

<t>r  [« 

4>n  [s 

4>t  C6 

1  SGI  [  7 

Irgi  [  8 
IN-Im  [ 9 
voi[io 
RESET  f  11 


20 J  Vdd 

19  ]  </>3~ 

18  ]  <£r3 

17  ]  tf>3+ 

16 ]  ISGR 
15  ]  IrGR 

14]  IN- Re 
13  ]  V0R 
12  ]  V0G 


Figure  1  -  Pin  Designations  for  R5601 


ISGR  IRGR 


OUTPUTS  TO 
DIFFERENTIAL 
AMPLIFIERS  I 


ISGI  IRGI 


Figure  2  -  Simplified  Block  Diagram  of  R5601 


^RGR, 
^  RGI 


1  SGR, 
I  SGI 


Figure  3  -  Clock  waveforms  and  timing  diagram. 


NOTES' 

1  component  values  are  TYPICAL. 

Z  "C"  IS  SAME  WAVEFORM  BUT  SEPARATE 
DRIVER  EG  <f>x  SAME  WAVEFORM  AS  <f>y 

3  IRGR  •  INPUT  RECV  GATE /REAL 
I  RGI  ■  INPUT  RECV  GATE/IMAG 
I  SGR  ■  INPUT  SAMPLE  GATE /REAL 
I  SGI  •  HPUT  SAMPLE  GATE/IMAG 
4.  (  I  INDICATES  TYPICAL  VOLTAGE  LEVEL 


(try) 


Figure  4  -  Connection  Arrangement  for  R5601. 
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